Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_000913 2,071,498 G→A 10.2% noncoding (182/195 nt) isrC → small RNA IsrC

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0009132,071,4980GA10.2% 282.2 / 15.9 98noncoding (182/195 nt)isrCsmall RNA IsrC
Reads supporting (aligned to +/- strand):  ref base G (39/49);  new base A (4/6);  total (43/55)
Fisher's exact test for biased strand distribution p-value = 1.00e+00
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 8.87e-01

CGATAAGCTAATAATAACCTTTGTCAGTAACATGCACAGATACGTACAGAAAGACATTCAGGGAACAACAGAACCACAATTCAGAAACTCCCACAGCCGGACCTCCGGCACTGTAACCCTTTACCTGCCGGTATCCACGTTTGTGGGTACCGGCTTTTTTATTCACCCTCAATCTAAGGAAAAGCTGATGAAACGACATCTGAATACCTGCTACAGGCTGGTATGGAATCACATGACGGGCGCTTTCGTGGTTGCCTCCGAACTGGCCCGCGCACGGGGTAAACGTGGCGGTG  >  NC_000913/2071352‑2071644
                                                                                                                                                  |                                                                                                                                                  
cgatAAGCTAATAATAACCTTTGTCAGTAACATGCACAGATACGTACAGAAAGACATTCAGGGAACAACAGAACCACAATTCAGAAACTCCCACAGCCGGACCTCCGGCACTGTAAGCCTTTACCTGCCGGTATCCACGTTTGTGGGTAc                                                                                                                                                 <  2:951540/150‑1 (MQ=255)
 gatAAGCTAATAATAACCTTTGTCAGTAACATGCACAGATACGTACAGAAAGACATTCAGGGAACAACAGAACCACAATTCAGAAACTCCCACAGCCGGACCTCCGGCACTGTAACCCTTTACCTGCCGGTATCCACGTTTGTGGGTAcc                                                                                                                                                <  1:125368/150‑1 (MQ=255)
   tAAGCTAATAATAACCTTTGTCAGTAACATGCACAGATACGTACAGAAAGACATTCAGGGAACAACAGAACCACAATTCAGAAACTCCCACAGCCGGACCTCCGGCACTGTAACCCTTTACCTGCCGGTACCCACAAACGTGGATACCgg                                                                                                                                              <  2:398940/150‑1 (MQ=255)
   tAAGCTAATAATAACCTTTGTCAGTAACATGCACAGATACGTACAGAAAGACATTCAGGGAACAACAGAACCACAATTCAGAAACTCCCACAGCCGGACCTCCGGCACTGTAACCCTTTACCTGCCGGTACCCACAAACGTGGATACCgg                                                                                                                                              <  1:1236015/150‑1 (MQ=255)
     aGCTAATAATAACCTTTGTCAGTAACATGCACAGATACGTACAGAAAGACATTCAGGGAACAACAGAACCACAATTCAGAAACTCCCACAGCCGGACCTCCGGCACTGTAACCCTTTACCTGCCGGTATCCACGTTTGTggg                                                                                                                                                    >  2:1390969/1‑142 (MQ=255)
     aGCTAATAATAACCTTTGTCAGTAACATGCACAGATACGTACAGAAAGACATTCAGGGAACAACAGAACCACAATTCAGAAACTCCCACAGCCGGACCTCCGGCACTGTAACCCTTTACCTGCCGGTATCCACGTTTGTggg                                                                                                                                                    <  1:1390969/142‑1 (MQ=255)
     aGCTAATAATAACCTTTGTCAGTAACATGCACAGATACGTACAGAAAGACATTCAGGGAACAACAGAACCACAATTCAGAAACTCCCACAGCCGGACCTCCGGCACTGTAACCCTTTACCTGCCGGTACCCACAAACGTGGATACCGGca                                                                                                                                            <  1:585683/150‑2 (MQ=255)
      gCTAATAATAACCTTTGTCAGTAACATGCACAGATACGTACAGAAAGACATTCAGGGAACAACAGAACCACAATTCAGAAACTCCCACAGCCGGACCTCCGGCACTGTAACCCTTTACCTGCCGGTACCCACAAACGTGGATACCGGcag                                                                                                                                           <  1:1297104/150‑3 (MQ=255)
      gCTAATAATAAACTTTGTCAGTAACATGCACAGATACGTACAGAAAGACATTAAGGGAAAAACAGAAACACAATTCAGAAACTCCCACAGCCGGACCTCCGGCACTGTAACCCTTTACCTGCCGGTATCCACGTTTGTGGGTACCGGCtt                                                                                                                                           >  2:207054/1‑150 (MQ=255)
         aataataaCCTTTGTCAGTAACATGCACAGATACGTACAGAAAGACATTCAGGGAACAACAGAACCACAATTCAGAAACTCCCACAGCCGGACCTCCGGCACTGTAACCCTTTACCTGCCGGTATCCACGTTTGTGGGTACCGGCttttt                                                                                                                                        <  1:1073210/150‑1 (MQ=255)
           taataaCCTTTGTCAGTAACATGCACAGATACGTACAGAAAGACATTCAGGGAACAACAGAACCACAATTCAGAAACTCCCACAGCCGGACCTCCGGCACTGTAACCCTTTACCTGCCGGTATCCACGTTTGTGGGTACCGGCTTTTTTa                                                                                                                                      <  1:483738/150‑1 (MQ=255)
           taataaCCTTTGTCAGTAACATGCACAGATACGTACAGAAAGACATTCAGGGAACAACAGAACCACAATTCAGAAACTCCCACAGCCGGACCTCCGGCACTGTAACCCTTTACCTGCCGGTATCCACGTTTGTGGGTACCGGCTTTTTTa                                                                                                                                      <  1:1580735/150‑1 (MQ=255)
             ataaCCTTTGTCAGTAACATGCACAGATACGTACAGAAAGACATTCAGGGAACAACAGAACCACAATTCAGAAACTCCCACAGCCGGACCTCCGGCACTGTAACCCTTTACCTGCCGGTATCCACGTTTGTGGGTACCGGCTTTTTTAtt                                                                                                                                    <  1:529451/150‑1 (MQ=255)
             ataaCCTTTGTCAGTAACATGCACAGATACGTACAGAAAGACATTCAGGGAACAACAGAACCACAATTCAGAAACTCCCACAGCCGGACCTCCGGCACTGTAACCCTTTACCTGCCGGTATCCACGTTTGTGGGTACCGGCTTTTTTAtt                                                                                                                                    <  1:228490/150‑1 (MQ=255)
              taaCCTTTGTCAGTAACATGCACAGATACGTACAGAAAGACATTCAGGGAACAACAGAACCACAATTCAGAAACTCCCACAGCCGGACCTCCGGCACTGTAACCCTTTACCTGCCGGTATCCACGTTTGTGGGTACCGGCTTTTTTATTc                                                                                                                                   <  2:362623/150‑1 (MQ=255)
                aCCTTTGTCAGTAACATGCACAGATACGTACAGAAAGACATTCAGGGAACAACAGAACCACAATTCAGAAACTCCCACAGCCGGACCTCCGGCACTGTAACCCTTTACCTGCCGGTATCCACGTTTGTGGGTACCGGCTTTTTTATTCAc                                                                                                                                 <  2:522150/150‑1 (MQ=255)
                 ccTTTGTCAGTAACATGCACAGATACGTACAGAAAGACATTCAGGGAACAACAGAACCACAATTCAGAAACTCCCACAGCCGGACCTCCGGCACTGTAACCCTTTACCTGCCGGTATCCACGTTTGTGGGTACCGGCTTTTTTATTCAcc                                                                                                                                <  2:531171/150‑1 (MQ=255)
                  cTTTGTCAGTAACATGCACAGATACGTACAGAAAGACATTCAGGGAACAACAGAACCACAATTCAGAAACTCCCACAGCCGGACCTCCGGCACTGTAACCCTTTACCTGCCGGTATCCACGTTTGTGGGTACCGGCTTTTTTATTCAccc                                                                                                                               <  1:1495375/150‑1 (MQ=255)
                   tttGTCAGTAACATGCACAGATACGTACAGAAAGACATTCAGGGAACAACAGAACCACAATTCAGAAACTCCCACAGCCGGACCTCCGGCACTGTAACCCTTTACCTGCCGGTATCCACGTTTGTGGGTACCGGCTTTTTTATTCACCCt                                                                                                                              <  2:240044/150‑1 (MQ=255)
                    ttGTCAGTAACATGCACAGATACGTACAGAAAGACATTCAGGGAACAACAGAACCACAATTCAGAAACTCCCACAGCCGGACCTCCGGCACTGTAACCCTTTACCTGCCGGTATCCACGTTTGTGGGTACCGGCTTTTTTATTCACCCTc                                                                                                                             <  1:461789/150‑1 (MQ=255)
                        cAGTAACATGCACAGATACGTACAGAAAGACATTAAGGGAACAACAGAACCACAATTCAGAAACTCCCACAGCCGGACCTCCGGCACTGTAACCCTTTACCTGCCGGTATCCACGTTTGTGGGTACCGGCTTTTTTATTCACCCTCAATc                                                                                                                         >  2:1494087/1‑150 (MQ=255)
                         aGTAACATGCACAGATACGTACAGAAAGACATTCAGGGAACAACAGAACCACAATTCAGAAACTCCCACAGCCGGACCTCCGGCACTGTAACCCTTTACCTGCCGGTATCCACGTTTGTGGGTACCGGCTTTTTTATTCACCCTCAATCt                                                                                                                        >  1:1496999/1‑150 (MQ=255)
                          gTAACATGCACAGATACGTACAGAAAGACATTCAGGGAACAACAGAACCACAATTCAGAAACTCCCACAGCCGGACCTCCGGCACTGTAACCCTTTACCTGCCGGTATCCACGTTTGTGGGTACCGGCTTTTTTATTCACCCTCAATCTa                                                                                                                       <  2:923071/150‑1 (MQ=255)
                            aaCATGCACAGATACGTACAGAAAGACATTCAGGGAACAACAGAACCACAATTCAGAAACTCCCACAGCCGGACCTCCGGCACTGTAACCCTTTACCTGCCGGTATCCACGTTTGTGGGTACCGGCTTTTTTATTCACACTCAATCTAAg                                                                                                                     >  2:318293/1‑150 (MQ=255)
                                tGCACAGATACGTACAGAAAGACATTCAGGGAACAACAGAACCACAATTCAGAAACTCCCACAGCCGGACCTCCGGCACTGTAACCCTTTACCTGCCGGTATCCACGTTTGTGGGTACCGGCTTTTTTATTCACCCTCAATCTAAGGaaa                                                                                                                 >  2:521185/1‑150 (MQ=255)
                                    caGATACGTACAGAAAGACATTCAGCGAACAACAGAACCACAATTCAGAAACTCCCACAGCCGGACCTCCGGCACTGTAACCCTTTACCTGCCGGTATCCACGTTTGTGGGTACCGGCTTTTTTATTCACCCTCAATCTAAGGAAAAGCt                                                                                                             >  2:1101263/1‑150 (MQ=255)
                                     aGATACGTACAGAAAGACATTCAGGGAACAACAGAACCACAATTCAGAAACTCCCACAGCCGGACCTCCGGCACTGTAACCCTTTACCTGCCGGTATCCACGTTTGTGGGTACCGGCTTTTTTATTCACCCTCAATCTAAGGAAAAGCtg                                                                                                            <  2:1008244/150‑1 (MQ=255)
                                       aTACGTACAGAAAGACATTCAGGGAACAACAGAACCACAATTCAGAAACTCCCACAGCCGGACCTCCGGCACTGTAACCCTTTACCTGCCGGTATCCACGTTTGTGGGTACCGGCTTTTTTATTCACCCTCAATCTAAGGAAAAGCtgat                                                                                                          <  2:1527984/150‑1 (MQ=255)
                                       aTACGTACAGAAAGACATTCAGGGAACAACAGAACCACAATTCAGAAACTCCCACAGCCGGACCTCCGGCACTGTAACCCTTTACCTGCCGGTATCCACGTTTGTGGGTACCGGCTTTTTTATTCACCCTCAATCTAAGGAAAAGCtgat                                                                                                          <  1:786010/150‑1 (MQ=255)
                                        tACGTACAGAAAGACGTTCAGGGAACAACAGAACCACAATTCAGAAACTCCCACAGCCGGACCTCCGGCACTGTAACCCTTTACCTGCCGGTATCCACGTTTGTGGGTACCGGCTTTTTTATTCACCCTCAATCTAAGGAAAAGCtgatg                                                                                                         <  2:1360219/150‑1 (MQ=255)
                                        tACGTACAGAAAGACATTCAGGGAACAACAGAACCACAATTCAGAAACTCCCACAGCCGGACCTCCGGCACTGTAACCCTTTACCTGCCGGTATCCACGTTTGTGGGTACCGGCTTTTTTATTCACCCTCAATCTAAGGAAAAGCtgatg                                                                                                         <  2:865466/150‑1 (MQ=255)
                                        tACGTACAGAAAGACATTCAGGGAACAACAGAACCACAATTCAGAAACTCCCACAGCCGGACCTCCGGCACTGTAACCCTTTACCTGCCGGTATCCACGTTTGTGGGTACCGGCTTTTTTATTCACCCTCAATCTAAGGAAAAGCtgatg                                                                                                         <  2:1060915/150‑1 (MQ=255)
                                          cGTACAGAAAGACATTCAGGGAACAACAGAACCACAATTCAGAAACTCCCACAGCCGCACCTCCGGCACAGTAACCCTTTACCTGCCGGTATCCACGTTTGTGGGTACCGGCTTTTTTATTCACCCTCAATCTAAGGAAAAGCTGATGaa                                                                                                       >  1:838474/1‑150 (MQ=255)
                                             aCAGAAAGACATTCAGGGAACAACAGAACCACAATTCAGAAACTCCCACAGCCGGACCTCCGGCACTGTAACCCTTTACCTGCCGGTATCCACGTTTGTGGGTACCGGCTTTTTTATTCACCCTCAATCTAAGGAAAAGCTGATGAAACg                                                                                                    >  1:427065/1‑150 (MQ=255)
                                              cAGAAAGACATTCAGGGAACAACAGAACCACAATTCAGAAACTCCCACAGCCGGACCTCCGGCACTGTAACCCTTTACCTGCCGGTATCCACGTTTGTGGGTACCGGCTTTTTTATTCACCCTCAATCTAAGGAAAAGCTGATGAAACGa                                                                                                   >  1:701607/1‑150 (MQ=255)
                                                gAAAGACATTCAGGGAACAACAGAACCACAATTCAGAAACTCCCACAGCCGGACCTCCGGCACTGTAACCCTTTACCTGCCGGTATCCACGTTTGTGGGTACCGGCTTTTTTATTCACCCTCAATCTAAGGAAAAGCTGATGAAACGACa                                                                                                 >  2:1285816/1‑150 (MQ=255)
                                                    gACATTCAGGGAACAACAGAACCACAATTCAGAAACTCCCACAGCCGGACCTCCGGCACTGTAACCCTTTACCTGCCGGTATCCACGTTTGTGGGTACCGGCTTTTTTATTCACCCTCAATCTAAGGAAAAGCTGATGAAACGACATCTg                                                                                             <  2:1523923/150‑1 (MQ=255)
                                                      cATTCAGGGAACAACAGAACCACAATTCAGAAACTCCCACAGCCGGACCTCCGGCACTGTAACCCTTTACCTGCCGGTATCCACGTTTGTGGGTACCGGCTTTTTTATTCACCCTCAATCTAAGGAAAAGCTGATGAAACGACATCTGaa                                                                                           >  1:892509/1‑150 (MQ=255)
                                                       aTTCAGGGAACAACAGAACCACAATTCAGAAACTCCCACAGCCGGACCTCCGGCACTGTAACCCTTTACCTGCCGGTATCCACGTTTGTGGGTACCGGCTTTTTTATTCACCCTCAATCTAAGGAAAAGCTGATGAAACGACATCTGAAt                                                                                          >  1:321361/1‑150 (MQ=255)
                                                       aTTCAGGGAACAACAGAACCACAATTCAGAAACTCCCACAGCCGGACCTCCGGCACTGTAACCCTTTACCTGCCGGTATCCACGTTTGTGGGTACCGGCTTTTTTATTCACCCTCAATCTAAGGAAAAGCTGATGAAACGACATCTGAAt                                                                                          >  1:1040917/1‑150 (MQ=255)
                                                           aGGGAACAACAGAACGACAATTCATAAACTCCCACAGCCGGACCTCCGGCACTGTAACGCTTTACCTGCCGGTATCCACGTTTGTGGGTACCGGCTTTTTTATTCACCCTCAATCTAAGGAAAAGCTGATGATACGACATCTGAATAcca                                                                                      >  2:1326904/1‑149 (MQ=255)
                                                           aGGGAACAACAGAACCACAATTCAGAAACTCCCACAGCCGGACCTCCGGCACTGTAACCCTTTACCTGCCGGTATCCACGTTTGTGGGTACCGGCTTTTTTATTCACCCTCAATCTAAGGAAAAGCTGATGAAACGACATCTGAATACCt                                                                                      <  2:410309/150‑1 (MQ=255)
                                                           aGGGAACAACAGAACCACAATTCAGAAACTCCCACAGCCGGACCTCCGGCACTGTAACCCTTTACCTGCCGGTATCCACGTTTGTGGGTACCGGCTTTTTTATTCACCCTCAATCTAAGGAAAAGCTGATGAAACGACATCTGAATACCt                                                                                      <  1:895331/150‑1 (MQ=255)
                                                              gAACAACATAACCACAATTCAGAAACTCCCACAGCCGGACCTCCGGCACTGTAACCCTTTACCTGCCGGTATCCACGTTTGTGGGTACCGGCTTTTTTATTCACCCTCAATCTAAGGAAAAGCTGATGAAACGACATCTGAATACCTGCt                                                                                   <  1:760013/150‑1 (MQ=255)
                                                                  aacaGAACCACAATTCAGAAACTCCCACAGCCGGACCTCCGGCACTGTAACCCTTTACCTGCCGGTATCCACGTTTGTGGGTACCGGCTTTTTTATTCACCCTCAATCTAAGGAAAAGCTGATGAAACGACATCTGAATACCTGCTACAg                                                                               >  2:243859/1‑150 (MQ=255)
                                                                  aacaGAACCACAATTCAGAAACTCCCACAGCCGGACCTCCGGCACTGTAACCCTTTACCTGCCGGTATCCACGTTTGTGGGTACCGGCTTTTTTATTCACCCTCAATCTAAGGAAAAGCTGATGAAACGACATCTGAATACCTGCTACAg                                                                               >  2:1391833/1‑150 (MQ=255)
                                                                    caGAACCACAATTCAGAAACTCCCACAGCCGGACCTCCGGCACTGTAACCCTTTACCTGCCGGTATCCACGTTTGTGGGTACCGGCTTTTTTATTCACCCTCAATCTAAGGAAAAGCTGATGAAACGACATCTGAATACCTGCTACAGGc                                                                             >  2:1418396/1‑150 (MQ=255)
                                                                     aGAACCACAATTCAGAAACTCCCACAGCCGGACCTCCGGCACTGTAACCCTTTACCTGCCGGTATCCACGTTTGTGGGTACCGGCTTTTTTATTCACCCTCAATCTAAGGAAAAGCTGATGAAACGACATCTGAATACCTGCTACAGGCt                                                                            >  1:172843/1‑150 (MQ=255)
                                                                      gAACCACAATTCAGAAACTCCCACAGCCGGACCTCCGGCACTGTAACCCTTTACCTGCCGGTATCCACGTTTGTGGGTACCGGCTTTTTTATTCACCCTCAATCTAAGGAAAAGCTGATGAAACGACATCTGAATACCTGCTACAGGCTg                                                                           >  1:494938/1‑150 (MQ=255)
                                                                       aaCCACAATTCAGAAACTCCCACAGCCGGACCTCCGGCACTGTAACCCTTTACCTGCCGGTATCCACGTTTGTGGGTACCGGCTTTTTTATTCACCCTCAATCTAAGGAAAAGCTGATGAAACGACATCTGAATACCTGCTACAGGCTgg                                                                          <  2:1533261/150‑1 (MQ=255)
                                                                       aaCCACAATTCAGAAACTCCCACAGCCGGACCTCCGGCACTGTAACCCTTTACCTGCCGGTATCCACGTTTGTGGGTACCGGCTTTTTTATTCACCCTCAATCTAAGGAAAAGCTGATGAAACGACATCTGAATACCTGCTACAGGCTgg                                                                          >  1:463964/1‑150 (MQ=255)
                                                                       aaCCACAATTCAGAAACTCCCACAGCCGGACCTCCGGCACTGTAACCCTTTACCTGCCGGTATCCACGTTTGTGGGTACCGGCTTTTTTATTCACCCTCAATCTAAGGAAAAGCTGATGAAACGACATCTGAATACCTGCTACAGGCTgg                                                                          <  1:216105/150‑1 (MQ=255)
                                                                        aCCACAATTCAGAAACTCCCACAGCCGGACCTCCGGCACTGTAACCCTTTACCTGCCGGTATCCACGTTTGTGGGTACCGGCTTTTTTATTCACCCTCAATCTAAGGAAAAGCTGATGAAACGACATCTGAATACCTGCTACAGGCTGGt                                                                         <  2:126673/150‑1 (MQ=255)
                                                                          cacaATTCAGAAACTCCCACAGCCGGACCTCCGGCACTGTAACCCTTTACCTGCCGGTATCCACGTTTGTGGGTACCGGCTTTTTTATTCACCCTCAATCTAAGGAAAAGCTGATGAAACGACATCTGAATACCTGCTACAGGCTGGTAt                                                                       >  2:36654/1‑150 (MQ=255)
                                                                           acaATTCAGAAACTCCCACAGCCGGACCTCCGGCACTGTAACCCTTTACCTGCCGGTATCCACGTTTGTGGGTACCGGCTTTTTTATTCACCCTCAATCTAAGGAAAAGCTGATGAAACGACATCTGAATACCTGCTACAGGCTGGTATg                                                                      >  2:1038082/1‑150 (MQ=255)
                                                                           acaATTCAGAAACTCCCACAGCCGGACCTCCGGCACTGTAACCCTTTACCTGCCGGTATCCACGTTTGTGGGTACCGGCTTTTTTATTCACCCTCAATCTAAGGAAAAGCTGATGAAACGACATCTGAATACCTGCTACAGGCTGGTATg                                                                      <  1:1101263/150‑1 (MQ=255)
                                                                            caATTCAGAAACTCCCACAGCCGGACCTCCGGCACTGTAACCCTTTACCTGCCGGTATCCACGTTTGTGGGTACCGGCTTTTTTATTCACCCTCAATCTAAGGAAAAGCTGATGAAACGACATCTGAATACCTGCTACAGGCTGGTATgg                                                                     >  2:793619/1‑150 (MQ=255)
                                                                                 cAGAAACTCCCACAGCCGGACCTCCGGCACTGTAACCCTTTACCTGCCGGTATCCACGTTTGTGGGTACCGGCTTTTTTATTCACCCTCAATCTAAGGAAAAGCTGATGAAACGACATCTGAATACCTGCTACAGGCTGGTATGGAATca                                                                <  1:901577/150‑1 (MQ=255)
                                                                                 cAGAAACTCCCACAGCCGGACCTCCGGCACTGTAACCCTTTACCTGCCGGTATCCACGTTTGTGGGTACCGGCTTTTTTATTCACCCTCAATCTAAGGAAAAGCTGATGAAACGACATCTGAATACCTGCTACAGGCTGGTATGGAATca                                                                <  1:1346070/150‑1 (MQ=255)
                                                                                 cAGAAACTCCCACAGCCGGACCTCCGGCACTGTAACCCTTTACCTGCCGGTATCCACGTTTGTGGGTACCGGCTTTTTTATTCACCCTCAATCTAAGGAAAAGCTGATGAAACGACATCTGAATACCTGCTACAGGCTGGTATGGAATca                                                                >  2:453219/1‑150 (MQ=255)
                                                                                   gAAACTCCCACAGCCGGACCTCCGGCACTGTAACCCTTTACCTGCCGGTATCCACGTTTGTGGGTACCGGCTTTTTTATTCACCCTCAATCTAAGGAAAAGCTGATGAAACGACATCTGAATACCTGCTACAGGCTGGTATGGAATcaca                                                              >  1:380007/1‑150 (MQ=255)
                                                                                   gAAACTCCCACAGCCGGACCTCCGGCACTGTAACCCTTTACCTGCCGGTATCCACGTTTGTGGGTACCGGCTTTTTTATTCACCCTCAATCTAAGGAAAAGCTGATGAAACGACATCTGAATACCTGCTACAGGCTGGTATGGAATcaca                                                              <  1:1418396/150‑1 (MQ=255)
                                                                                     aaCTCCCACAGCCGGACCTCCGGCACTGTAACCCTTTACCTGCCGGTATCCACGTTTGTGGGTACCGGCTTTTTTATTCACCCTCAATCTAAGGAAAAGCTGATGAAACGACATCTGAATACCTGCTACAGGCTGGTTTGGAATCACATg                                                            >  1:1162739/1‑150 (MQ=255)
                                                                                       cTCCCACAGCCGGACCTCCGGCACTGTAACCCTTTACCTGCCGGTATCCACGTTTGTGGGTACCGGCTTTTTTATTCACCCTCAATCTAAGGAAAAGCTGATGAAACGACATCTGAATACCTGCTACAGGCTGGTATGGAATCACATGAc                                                          >  1:705759/1‑150 (MQ=255)
                                                                                        tCCCACAGCCGGACCTCCGGCACTGTAACCCTTTACCTGCCGGTATCCACGTTTGTGGGTACCGGCTTTTTTATTCACCCTCAATCTAAGGAAAAGCTGATGAAACGACATCTGAATACCTGCTACAGGCTGGTATGGAATCACATGACg                                                         >  1:488706/1‑150 (MQ=255)
                                                                                           cagtgCCGGAGGTCCGGCACTGTAACCCTTTACCTGCCGGTATCCACGTTTGTGGGTACCGGCTTTTTTATTCACCCTCAATCTAAGGAAAAGCTGATGAAACGACATCTGAATACCTGCTACAGGCTGGTATGGAATCACATGACGGgc                                                      >  2:1535622/5‑150 (MQ=255)
                                                                                             caGCCGGACCTCCGGCACTGTAACCCTTTACCTGCCGGTATCCACGTTTGTGGGTACCGGCTTTTTTATTCACCCTCAATCTAAGGAAAAGCTGATGAAACGACATCTGAATACCTGCTACAGGCTGGTATGGAATCACATGACGGgcgc                                                    >  2:341861/1‑150 (MQ=255)
                                                                                               gCCGGACCTCCGGCACTGTAACCCTTTACCTGCCGGTATCCACGTTTGTGGGTACCGGCTTTTTTATTCACCCTCAATCTAAGGAAAAGCTGATGAAACGACATCTGAATACCTGCTACAGGCTGGTATGGAATCACATGACGGGCGCtt                                                  >  2:109285/1‑150 (MQ=255)
                                                                                                   gACCTCCGGCACTGTAACCCTTTACCTGCCGGTATCCACGTTTGTGGGTACCGGCTTTTTTATTCACCCTCAATCTAAGGAAAAGCTGATGAAACGACATCTGAATACCTGCTACAGGCTGGTATGGAATCACATGACGGGCGCTTTCGt                                              <  1:1495933/150‑1 (MQ=255)
                                                                                                     ccTCCGGCACTGTAACCCTTTACCTGCCGGTATCCACGTTTGTGGGTACCGGCTTTTTTATTCACCCTCAATCTAAGGAAAAGCTGATGAAACGACATCTGAATACCTGCTACAGGCTGGTATGGAATCACATGACGGGCGCTTTCGTgg                                            <  1:168159/150‑1 (MQ=255)
                                                                                                     ccTCCGGCACTGTAACCCTTTACCTGCCGGTATCCACGTTTGTGGGTACCGGCTTTTTTATTCACCCTCAATCTAAGGAAAAGCTGATGAAACGACATCTGAATACCTGCTACAGGCTGGTATGGAATCACATGACGGGCGCTTTCGTgg                                            >  1:347069/1‑150 (MQ=255)
                                                                                                      cTCCGGCACTGTAACCCTTTACCTGCCGGTATCCACGTTTGTGGGTACCGGCTTTTTTATTCACCCTCAATCTAAGGAAAAGCTGATGAAACGACATCTGAATACCTGCTACAGGCTGGTATGGAATCACATGACGGGCGCTTTCGTGGt                                           <  1:1612941/150‑1 (MQ=255)
                                                                                                       tCCGGCACTGTAACCCTTTACCTGCCGGTATCCACGTTTGTGGGTACCGGCTTTTTTATTCACCCTCAATCTAAGGAAAAGCTGATGAAACGACATCTGAATACCTGCTACAGGCTGGTATGGAATCACATGACGGGCGCTTTCGTGGtt                                          <  1:207054/150‑1 (MQ=255)
                                                                                                        ccGGCACTGTAACCCTTTACCTGCCGGTATCCACGTTTGTGGGTACCGGCTTTTTTATTCACCCTCAATCTAAGGAAAAGCTGATGAAACGACATCTGAATACCTGCTACAGGCTGGTATGGAATCACATGACGGGCGCTTTCGTGGTTg                                         >  1:801275/1‑150 (MQ=255)
                                                                                                         cGGCACTGTAACCCTTTACCTGCCGGTATCCACGTTTGTGGGTACCGGCTTTTTTATTCACCCTCAATCTAAGGAAAAGCTGATGAAACGACATCTGAATACCTGCTACAGGCTGGTATGGAATCACATGACGGGCGCTTTCGTGGTTGc                                        >  1:1502991/1‑150 (MQ=255)
                                                                                                             aCTGTAACCCTTTACCTGCCGGTATCCACGTTTGTGGGTACCGGCTTTTTTATTCACCCTCAATCTAAGGAAAAGCTGATGAAACGACATCTGAATACCTGCTACAGGCTGGTATGGAATCACATGACGGGCGCTTTCGTGGTTGCCTcc                                    <  2:1108081/150‑1 (MQ=255)
                                                                                                                 tAACCCTTTACCTGCCGGTATCCACGTTTGTGGGTACCGGCTTTTTTATTCACCCTCAATCTAAGGAAAAGCTGATGAAACGACATCTGAATACCTGCTACAGGCTGGTATGGAATCACATGACGGGCGCTTTCGTGGTTGCCTCCGAAc                                <  1:1317809/150‑1 (MQ=255)
                                                                                                                     ccTTTACCTGCCGGTATCCACGTTTGTGGGTACCGGCTTTTTTATTCACCCTCAATCTAAGGAAAAGCTGATGAAACGACATCTGAATACCTGCTACAGGCTGGTATGGAATCACATGACGGGCGCTTTCGTGGTTGCCTCCGAACTGGc                            >  1:1225985/1‑150 (MQ=255)
                                                                                                                        ttACCTGCCGGTATCCACGTTTGTGGGTACCGGCTTTTTTATTCACCCTCAATCTAAGGAAAAGCTGATGAAACGACATCTGAATACCTGCTACAGGCTGGTATGGAATCACATGACGGGCGCTTTCGTGGTTGCCTCCGAACTGGCCcg                         <  2:347069/150‑1 (MQ=255)
                                                                                                                         tACCTGCCGGTATCCACGTTTGTGGGTACCGGCTTTTTTATTCACCCTCAATCTAAGGAAAAGCTGATGAAACGACATCTGAATACCTGCTACAGGCTGGTATGGAATCACATGACGGGCGCTTTCGTGGTTGCCTCCGAACTGGCCcgc                        >  1:973175/1‑150 (MQ=255)
                                                                                                                           aaagCCGGTACCCACAAACGTGGATACCGGCTTTTTTATTCACCCTCAATCTAAGGAAAAGCTGATGAAACGACATCTGAATACCTGCTACAGGCTGGTATGGAATCACATGACGGGCGCTTTCGTGGTTGCCTCCGAACTGGCCCgcgc                      >  2:1231363/4‑150 (MQ=255)
                                                                                                                            aagCCGGTACCCACAAACGTGGATACCGGCTTTTTTATTCACCCTCAATCTAAGGAAAAGCTGATGAAACGACATCTGAATACCTGCTACAGGCTGGTATGGAATCACATGACGGGCGCTTTCGTGGTTGCCTCCGAACTGGCCCgcgcc                     >  1:1353220/3‑149 (MQ=255)
                                                                                                                             agCCGGTACCCACAAACGTGGATACCGGCTTTTTTATTCACCCTCAATCTAAGGAAAAGCTGATGAACCGACATCTGAATACCTGCTACAGGCTGGTATGGAATCACATGACGGGCGCTTTCGTGGGTGCCTCCGAACTGGCCCgcgcgc                    >  1:719607/2‑148 (MQ=255)
                                                                                                                              gCCGGTATCCACGTTTGTGGGTACCGGCTTTTTTATTCACCCTCAATCTAAGGAAAAGCTGATGAAACGACATCTGAATACCTGCTACAGGCTGGTATGGAATCACATGACGGGCGCTTTCGTGGTTGCCTCCGAACTGGCCCGCGCACg                   >  1:1022936/1‑150 (MQ=255)
                                                                                                                              gCCGGTACCCACAAACGTGGATACCGGCTTTTTTATTCACCCTCAATCTAAGGAAAAGCTGATGAAACGACATCTGAATACCTGCTACAGGCTGGTATGGAATCACATGACGGGCGCTTTCGTGGTTg                                         <  1:586864/128‑1 (MQ=255)
                                                                                                                              gCCGGTACCCACAAACGTGGATACCGGCTTTTTTATTCACCCTCAATCTAAGGAAAAGCTGATGAAACGACATCTGAATACCTGCTACAGGCTGGTATGGAATCACATGACGGGCGCTTTCGTGGTTg                                         >  2:586864/1‑128 (MQ=255)
                                                                                                                               ccGGTATCCACGTTTGTGGGTACCGGCTTTTTTATTCACCCTCAATCTAAGGAAAAGCTGATGAAACGACATCTGAATACCTGCTACAGGCTGGTATGGAATCACATGACGGGCGCTTTCGTGGTTGCCTCCGAACTGGCCCGCGCACgg                  >  1:1458921/1‑150 (MQ=255)
                                                                                                                                cGGTATCCACGTTTGTGGGTACCGGCTTTTTTATTCACCCTCAATCTAAGGAAAAGCTGATGAAACGACATCTGAATACCTGCTACAGGCTGGTATGGAATCACATGACGGGCGCTTTCGTGGTTGCCTCCGAACTGGCCCGCGCACggg                 <  2:890488/150‑1 (MQ=255)
                                                                                                                                cGGTATCCACGTTTGTGGGTACCGGCTTTTTTATTCACCCTCAATCTAAGGAAAAGCTGATGAAACGACATCTGAATACCTGCTACAGGCTGGTATGGAATCACATGACGGGCGCTTTCGTGGTTGCCTCCGAACTGGCCCGCGCACggg                 <  1:109285/150‑1 (MQ=255)
                                                                                                                                 ggTATCCACGTTTGTGGGTACCGGCTTTTTTATTCACCCTCAATCTAAGGAAAAGCTGATGAAACGACATCTGAATACCTGCTACAGGCTGGTATGGAATCACATGACGGGCGCTTTCGTGGTTGCCTCCGAACTGGCCCGCGCACgggg                <  1:867788/150‑1 (MQ=255)
                                                                                                                                 ggTATCCACGTTTGTGGGTACCGGCTTTTTTATTCACCCTCAATCTAAGGAAAAGCTGATGAAACGACATCTGAATACCTGCTACAGGCTGGTATGGAATCACATGACGGGCGCTTTCGTGGTTGCCTCCGAACTGGCCCGCGCACgggg                <  1:1285816/150‑1 (MQ=255)
                                                                                                                                  gTATCCACGTTTGTGGGTACCGGCTTTTTTATTCACCCTCAATCTAAGGAAAAGCTGATGAAACGACATCTGAATACCTGCTACAGGCTGGTATGGAATCACATGACGGGCGCTTTCGTGGTTGCCTCCGAACTGGCCCGCGCACGGGGt               <  1:1038082/150‑1 (MQ=255)
                                                                                                                                  gTATCCACGTTTGTGGGTACCGGCTTTTTTATTCACCCTCAATCTAAGGAAAAGCTGATGAAACGACATCTGAATACCTGCTACAGGCTGGTATGGAATCACATGACGGGCGCTTTCGTGGTTGCCTCCGAACTGGCCCGCGCACGGGGt               >  1:1406649/1‑150 (MQ=255)
                                                                                                                                  gTATCCACGTTGTGGGTACCGGCTTTTTTATTCACCCTCAATCTAAGGAAAAGCTGATGAAACGACATCTGAATACCTGCTACAGGCTGGTATGGAATCACATGACGGGCGCTTTCGTGGTTGCCTCCGAACTGGCCCGCGCACGGGGTa              <  1:639188/150‑1 (MQ=255)
                                                                                                                                     tCCACGTTTGTGGGTACCGGCTTTTTTATTCACCCTCAATCTAAGGAAAAGCTGATGAAACGACATCTGAATACCTGCTACAGGCTGGTATGGAATCACATGACGGGCGCTTTCGTGTTTGCCTCCGAACTGGCCCGCGCACGGGGTaaa            <  2:939288/150‑1 (MQ=255)
                                                                                                                                         cGTTTGTGGGTACCGGCTTTTTTATTCACCCTCAATCTAAGGAAAAGCTGATGAAACGACATCTGAATACCTGCTACAGGCTGGTATGGAATCACATGACGGGCGCTTTCGTGGTTGCCTCCGAACTGGCCCGCGCACGGGGTAAACgtg        <  1:793619/150‑1 (MQ=255)
                                                                                                                                           tttGTGGGTACCGGCTTTTTTATTCACCCTCAATCTAAGGAAAAGCTGATGAAACGACATCTGAATACCTGCTACAGGCTGGTATGGAATCACATGACGGGCGCTTTCGTGGTTGCCTCCGAACTGGCCCGCGCACGGGGTAAACgtggc      <  2:120894/150‑1 (MQ=255)
                                                                                                                                              gtgGGTACCGGCTTTTTTATTCACCCTCAATCTAAGGAAAAGCTGATGAAACGACATCTGAATACCTGCTACAGGCTGGTATGGAATCACATGACGGGCGCTTTCGTGGTTGCCTCCGAACTGGCCCGCGCACGGGGTAAACgtggcggt   <  1:36654/150‑1 (MQ=255)
                                                                                                                                               tgGGTACCGGCTTTTTTATTCACCCTCAATCTAAGGAAAAGCTGATGAAACGACATCTGAATACCTGCTACAGGCTGGTATGGAATCACATGACGGGCGCTTTCGTGGTTGCCTCCGAACTGGCCCGCGCACGGGGTAAACgtggcggtg  <  2:494938/150‑1 (MQ=255)
                                                                                                                                               tgGATACCGGCTTTTTTATTCACCCTCAATCTAAGGAAAAGCTGATGAAACGACATCTGAATACCTGCTACAGGCTGGTATGGAATCACATGACGGGCGCTTTCGTGGTTGCCTCCGAACTGGCCCGCGCACGGGGTAAACgtggcggtg  <  1:1459846/150‑1 (MQ=255)
                                                                                                                                                  |                                                                                                                                                  
CGATAAGCTAATAATAACCTTTGTCAGTAACATGCACAGATACGTACAGAAAGACATTCAGGGAACAACAGAACCACAATTCAGAAACTCCCACAGCCGGACCTCCGGCACTGTAACCCTTTACCTGCCGGTATCCACGTTTGTGGGTACCGGCTTTTTTATTCACCCTCAATCTAAGGAAAAGCTGATGAAACGACATCTGAATACCTGCTACAGGCTGGTATGGAATCACATGACGGGCGCTTTCGTGGTTGCCTCCGAACTGGCCCGCGCACGGGGTAAACGTGGCGGTG  >  NC_000913/2071352‑2071644

Alignment Legend
Aligned base mismatch/match (shaded by quality score): ATCG/ATCG < 3 ≤ ATCG/ATCG < 14 ≤ ATCG/ATCG < 32 ≤ ATCG/ATCG < 36 ≤ ATCG/ATCG
Unaligned base: atcg    Masked matching base: atcg    Alignment gap:     Deleted base: