Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_000913 2,518,197 C→G 6.5% noncoding (35/76 nt) alaW ← tRNA‑Ala

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0009132,518,1970CG 6.5% 257.2 / 7.1 77noncoding (35/76 nt)alaWtRNA‑Ala
Reads supporting (aligned to +/- strand):  ref base C (46/26);  new base G (3/2);  total (49/28)
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.04e-01

CTAAGCGGGATCGAACCGCTGACCTCTTGCATGCCATGCAAGCGCTCTCCCAGCTGAGCTATAGCCCCACGATGCGTTTACGTACCAAGTTTGCTGGGTGCAAAATTTGGTGGAGCTAAGCGGGATCGAACCGCTGACCTCTTGCATGCCATGCAAGCGCTCTCCCAGCTGAGCTATAGCCCCGTCACGTAAAGCTTGTCGAGTTGACGGGCGGCA‑TCATATGAATTCCGCCCGAATGTGTCAACGGCAAATTGCAACGTGTAGTTTCAATCGCTGAAAAATCAGGCAAATGAACAAT  >  NC_000913/2518049‑2518346
                                                                                                                                                    |                                                                                                                                                      
cTAAGGGGGATCGAACCGCTGACCTCTTGCATGCCATGCTAGCGCTCTCCCAGCTGAGCTATAGCCCCACGATGCGTTTACGTACCAAGTTTGGTGGGTGCAAAATTTGGTGGAGCTAAGCGGGATCGAACCGCTGACCTCTTGcatgcc                                                                                                                                                       <  2:18244/150‑1 (MQ=255)
cTAAGCGGGATCGAACCGCTGACCTCTTGCATGCCATGCAAGCGCTCTCCCAGCTGAGCTATAGCCCCACGATGCGTTTACGTACCAAGTTTGCTGGGTGCAAAATTTGGTGGAGCTAAGCGGGATCGAACCGCTGACCTCTTGcatgcc                                                                                                                                                       >  1:230/1‑150 (MQ=255)
         aTCTAACCCCTGAGCTCTTCCCTGCCATGATAGCGCTCTCCCATCTGTGCTATATTCACACTATTAGTTTACGTACCATGTTTTCAGGGTTCTAAATTTGGTGGAGCTAAGCGGGATCGAACCGCTGACCTCTTGCATGGCATGCAAgag                                                                                                                                              <  2:909493/150‑3 (MQ=12)
          tCGAACCGCTGACCTCTTGCATGCCATGCAAGCGCTCTCCCAGCTGAGCTATAGCCCCACGATGCCTTTACGTACCAAGTTTGCTGGGTGCAAAATTTGGTGGAGCTAAGCGGGATCGAACCGCTGACCTCTTGCATGCCATGCAAgcgc                                                                                                                                             <  2:47006/150‑1 (MQ=255)
              aCCTCTGACCTCTTTCATGCCATGCAAGCGCTCTCCCATCTGAGCTATAGCCCCACGATGCGATTACGTACCAAGTTTGCTGGGTGCAAAATTTGGTGGAGCTAAGCGGGATCGAACCGCTGACCTCTTGCATGCCATGCAAGCGCtctc                                                                                                                                         <  2:611156/150‑1 (MQ=255)
               ccGCTGACCTCTTGCAGGCCCTGCTAGCGCTCTCCCAGCTGAGCTATACCCCCACGATGCGTTTACGTACCAAGTTTGCTGGGTGCAAAATTTGGTGGAGCTAAGCGGGATCGAACCGCTGACCTCTTGCATGCCATGCAAGCGCTCTcc                                                                                                                                        <  1:419098/150‑1 (MQ=255)
                        tctTGCATGCCATGCAAGCGCTCTCCCAGCTGAGCTATAGCCCCACGATGCGTTTACGTACCAAGTTTGCTGGGTGCAAAATTTGGTGGAGCTAAGCGGGATCGAACCGCTGACCTCTTGCATGCCATGCAAGCGCTCTCCCAGCTGAGc                                                                                                                               <  1:219050/150‑1 (MQ=255)
                                     gcAAGCGCTCTCCCAGCTGAGCTATAGCCCCACGATGCGTTTACGTACCAAGTTTGCTGGGTGCAAAATTTGGTGGAGCTAAGCGGGATCGAACCGCTGACCTCTTGCATGCCATGCAAgcg                                                                                                                                              <  1:1013484/122‑1 (MQ=255)
                                     gcAAGCGCTCTCCCAGCTGAGCTATAGCCCCACGATGCGTTTACGTACCAAGTTTGCTGGGTGCAAAATTTGGTGGAGCTAAGCGGGATCGAACCGCTGACCTCTTGCATGCCATGCAAgcg                                                                                                                                              >  2:1013484/1‑122 (MQ=255)
                                                  cAGCTGAGCTATAGCCCCACGATGCGTTTACGTACCAAGTTTGCTGGGTGCAAAATTTGGTGGAGCTAAGCGGGATCGAACCGCTGACCTCTTGCATGCCATGCAAGCGCTCTCCCAGCTGAGCTATAGCACCGTAACGTAAAGCTTGTc                                                                                                     >  2:59205/1‑150 (MQ=255)
                                                   atcTGAGCTATAGGCCCAGTATGCTTTTACGTTCCAAGTTTGCTGGGTGGAAAATTTGGTGGAGCTAAGCGGGATCGAACCGCTGACCTCTTGCATGCCATGCAAGCGCTCTCCCAGCTGAGCTATAGCCCCGTCACGTAAAGCTTGTCg                                                                                                    <  1:1017565/148‑1 (MQ=21)
                                                    gCTGAGCTATAGCCCCACGATGCGTTTACGTACCAAGTTTGCTGGGTGCAAAATTTGGTGGAGCTAAGCGGGATCGAACCGCTGACCTCTTGCATGCCATGCAAGCGCTCTCCCAGCTGAGCTATAGCCCCGTCACGTAAAGCTTGTCGa                                                                                                   >  2:285413/1‑150 (MQ=255)
                                                         gCTATAGCCCCACGATGCGTTTACGTACCAAGTTTGCTGGGTGCAAAATTTGGTGGAGCTAAGCGGGATCGAACCGCTGACCTCTTGCATGGCATGCAAg                                                                                                                                                >  1:909863/1‑100 (MQ=255)
                                                         gCTATAGCCCCACGATGCGTTTACGTACCAAGTTTGCTGGGTGCAAAATTTGGTGGAGCTAAGCGGGATCGAACCGCTGACCTCTTGCATGGCATGCAAg                                                                                                                                                <  2:909863/100‑1 (MQ=255)
                                                          cTATAGCCCCACGATGCGTTTACGTACCAAGTTTGCTGGGTGCAAAATTTGGTGGAGCTAAGCGGGATCGAACCGCTGACCTCTTGCATGCCATGCAAGCGCTCTCCCAGCTGAGCTATAGCCCCGTCACGTAAAGCTTGTCGAGTTGAc                                                                                             >  1:626602/1‑150 (MQ=255)
                                                               gCCCCACGATGCGTTTACGTACCAAGTTTGCTGGGTGCAAAATTTGGTGGAGCTAAGCGGGATCGAACCGCTGACCTCTTGCATGCCATGCAAGCGCTCTCCCAGCTGAGCTATAGCCCCGTCACGTAAAGCTTGTCGAGTTGACGggc                                                                                         <  2:1611664/149‑1 (MQ=255)
                                                               gCCCCACGATGCGTTTACGTACCAAGTTTGCTGGGTGCAAAATTTGGTGGAGCTAAGCGGGATCGAACCGCTGACCTCTTGCATGCCATGCAAGCGCTCTCCCAGCTGAGCTATAGCCCCGTCACGTAAAGCTTGTCGAGTTGACGggc                                                                                         >  1:1611664/1‑149 (MQ=255)
                                                                      gATGCGTTTACGTACCAAGTTTGCTGGGTGCAAAATTTGGTGGAGCTAAGCGGGATCGAACCGCTGACCTCTTGCATGCCATGCAAGCGCTCTCCCAGCTGAGCTATAGCCCCGTCACGTAAAGCTTGTCGAGTTGACGGGCGGca‑tcat                                                                                >  1:337635/1‑150 (MQ=255)
                                                                        tGCGTTTACGTACCAAGTTTGCTGGGTGCAAAATTTGGTGGAGCTAAGCGGGATCGAACCGCTGACCTCTTGCATGCCATGCAAGCGCTCTCCCAGCTGAGCTATAGCCCCGTCACGTAAAGCTTGTCGAGTTGACGGGCGGCA‑TCatat                                                                              >  1:134648/1‑150 (MQ=255)
                                                                        tGCGTTTACGTACCAAGTTTGCTGGGTGCAAAATTTGGTGGAGCTAAGCGGGATCGAACCGCTGACCTCTTGCATGCCATGCAAGCGCTCTCCCAGCTGAGCTATAGCCCCGTCACGTAAAGCTTGTCGAGTTGACGGGCGGCA‑TCatat                                                                              >  2:679068/1‑150 (MQ=255)
                                                                          cGTTTACGTACCAAGTTTGCTGGGTGCAAAATTTGGTGGAGCTAAGCGGGATCGAACCGCTGACCTCTTGCATGCCATGCAAGCGCTCTCCCAGCTGAGCTATAGCCCCGTCACGTAAAGCTTGTCGAGTTGACGGGCGGCA‑TCATATGa                                                                            >  2:79174/1‑150 (MQ=255)
                                                                           gTTTACGTACCAAGTTTGCTGGGTGCAAAATTTGGTGGAGCTAAGCGGGATCGAACCGCTGACCTCTTGCATGCCATGCAAGCGCTCTCCCAGCTGAGCTATAGCCCCGTCACGTAAAGCTTGTCGAGTTGACGGGCGGCA‑TCATATGaa                                                                           >  1:979809/1‑150 (MQ=255)
                                                                                cGTACCAAGTTTGCTGGGTGCAAAATTTGGTGGAGCTAAGCGGGATCGAACCGCTGACCTCTTGCATGCCATGCAAGCGCTCTCCCAGCTGAGCTATAGCCCCGTCACGTAAAGCTTGTCGAGTTGACGGGCGGCA‑TCATATGAATTCCg                                                                      >  1:568147/1‑150 (MQ=255)
                                                                                 gTACCAAGTTTGCTGGGTGCAAAATTTGGTGGAGCTAAGCGGGATCGAACCGCTGACCTCTTGCATGCCATGCAAGCGCTCTCCCAGCTGAGCTATAGCCCCGTCACGTAAAGCTTGTCGAGTTGACGGGCGGCA‑TCATATGAATTCCGc                                                                     >  1:1354298/1‑150 (MQ=255)
                                                                                   aCCAAGTTTGCTGGGTGCAAAATTTGGTGGAGCTAAGCGGGATCGAACCGCTGACCTCTTTCATGCCATGCAAGCGCTCTCCCAGCTGAGCTATAGCCCCGTCACGTAAAGCTTGTCGAGTTGACGGGCGGCA‑TCATATGAATTCCGccc                                                                   >  1:827524/1‑150 (MQ=255)
                                                                                   aCCAAGTTTGCTGGGTGCAAAATTTGGTGGAGCTAAGCGGGATCGAACCGCTGACCTCTTGCATGCCATGCAAGCGCTCTCCCAGCTGAGCTATAGCCCCGTCACGTAAAGCTTGTCGAGTTGACGGGCGGCA‑TCATATGAATTCCGccc                                                                   >  2:163794/1‑150 (MQ=255)
                                                                                   aCCAAGTTTGCTGGGTGCAAAATTTGGTGGAGCTAAGCGGGATCGAACCGCTGACCTCTTGCATGCCATGCAAGCGCTCTCCCAGCTGAGCTATAGCCCCGTCACGTAAAGCTTGTCGAGTTGACGGGCGGCA‑TCATATGAATTCCGccc                                                                   >  2:274622/1‑150 (MQ=255)
                                                                                      aaGTTTGCTGGGTGCAAAATTTGGTGGAGCTAAGCGGGATCGAACCGCTGACCTCTTGCATGCCATGCAAGCGCTCTCCCAGCTGAGCTATAGCCCCGTCACGTAAAGCTTGTCGAGTTGACGGGCGGCA‑TCatat                                                                              <  2:1536879/136‑1 (MQ=255)
                                                                                      aaGTTTGCTGGGTGCAAAATTTGGTGGAGCTAAGCGGGATCGAACCGCTGACCTCTTGCATGCCATGCAAGCGCTCTCCCAGCTGAGCTATAGCCCCGTCACGTAAAGCTTGTCGAGTTGACGGGCGGCA‑TCatat                                                                              >  1:1536879/1‑136 (MQ=255)
                                                                                      aaGTTTGCTGGGTGCAAAATTTGGTGGAGCTAAGCGGGATCGAACCGCTGACCTCTTGCATGCCATGCAAGCGCTCTCCCAGCTGAGCTATAGCCCCGTCACGTAAAGCTTGTCGAGTTGACGGGCGGCA‑TCATATGAATTCCGCCCGaa                                                                >  1:1396564/1‑150 (MQ=255)
                                                                                      aaGTTTGCTGGGTGCAAAATTTGGTGGAGCTAAGCGGGATCGAACCGCTGACCTCTTGCATGCCATGCAAGCGCTCTCCCAGCTGAGCTATAGCCCCGTCACGTAAAGCTTGTCGAGTTGACGGGCGGCA‑TCATATGAATTCCGCCCGaa                                                                <  1:610292/150‑1 (MQ=255)
                                                                                      aaGTTTGCTGGGTGCAAAATTTGGTGGAGCTAAGCGGGATCGAACCGCTGACCTCTTGCATGCCATGCAAGCGCTCTCCCAGCTGAGCTATAGCCCCGTCACGTAAAGCTTGTCGAGTTGACGGGCGGAATTCATATG‑ATGCCGCCCgtc                                                                <  1:51553/150‑3 (MQ=21)
                                                                                        gcaggCTGGGTGCAAAATTTGGTGGAGCTAAGCGGGATCGAACCGCTGACCTCTTGCATGCCATGCAAGCGCTCTCCCAGCTGAGCTATAGCCCCGTCACGTAAAGCTTGTCGAGTTGACGGGCGGCA‑TCATATGAATTCCGCCCGAAtg                                                              <  1:1213548/146‑1 (MQ=255)
                                                                                            gCTGGGTGCAAAATTTGGTGGAGCTAAGCGGGATCGAACCGCTGACCTCTTGCATGCCATGCAAGCGCTCTCCCAGCTGAGCTATAGCCCCGTCACGTAAAGCTTGTCGAGTTGACGGGCGGCA‑TCATATGAATTCCGCCCGAAtg                                                              >  2:281406/1‑146 (MQ=255)
                                                                                            gCTGGGTGCAAAATTTGGTGGAGCTAAGCGGGATCGAACCGCTGACCTCTTGCATGCCATGCAAGCGCTCTCCCAGCTGAGCTATAGCCCCGTCACGTAAAGCTTGTCGAGTTGACGGGCGGCA‑TCATATGAATTCCGCCCGAAtg                                                              <  1:281406/146‑1 (MQ=255)
                                                                                             cTGGGTGCAAAATTTGGTGGAGCTAAGCGGGATCGAACCGCTGACCTCTTGCATGCCATGCAAGCGCTCTCCCAGCTGAGCTATAGCCCCGTCACGTAAAGCTTGTCGAGTTGACGGGCGGCA‑TCATATGAATTCCGCCCGAATGTGTCa                                                         >  2:916125/1‑150 (MQ=255)
                                                                                               gggTGCAAAATTTGGTGGAGCTAAGCGGGATCGAACCGCTGACCTCTTGCATGCCATGCAAGCGCTCTCCCAGCTGAGCTATAGCCCCGTCACGTAAAGCTTGTCGAGTTGACGGGCGGCA‑TCATATGAATTCCGCCCGAATGTGTCAAc                                                       >  1:391645/1‑150 (MQ=255)
                                                                                                  tGCAAAATTTGGTGGAGCTAAGCGGGATCGAACCGCTGACCTCTTGCATGCCATGCAAGCGCTCTCCCAGCTGAGCTATAGCCCCGTCACGTAAAGCTTGTCGAGTTGACGGGCGGCA‑TCATATGAATTCCGCCCGAATGTGTCAACGGc                                                    >  2:151100/1‑150 (MQ=255)
                                                                                                   gCAAAATTTGGTGGAGCTAAGCGGGATCGAACCGCTGACCTCTTGCATGCCATGCAAGCGCTCTCCCAGCTGAGCTATAGCCCCGTCACGTAAAGCTTGTCGAGTTGACGGGCGGCA‑TCATATGAATTCCGCCCGAATGTGTCAACGGCa                                                   >  2:949878/1‑150 (MQ=255)
                                                                                                     aaaaTTTGGTGGAGCTAAGCGGGATCGAACCGCTGACCTCTTGCATGCCATGCAAGCGCTCTCCCAGCTGAGCTATAGCCCCGTCACGTAAAGCTTGTCGAGTTGACGGGCGGCA‑TCATATGAATTCCGCCCGAATGTGTCAACGGCaaa                                                 >  2:225082/1‑150 (MQ=255)
                                                                                                          ttGGTGGAGCTAAGCGGGATCGAACCGCTGACCTCTTGCATGCCATGCAAGCGCTCTCCCAGCTGAGCTATAGCCCCGTCACGTAAAGCTTGTCGAGTTGACGGGCGGCA‑TCATATGAATTCCGCCCGAATGTGTCAACGGCAAATTGCa                                            <  2:134648/150‑1 (MQ=255)
                                                                                                           tggtggAGCTAAGCGGGATCGAACCGCTGACCTCTTGCATGCCATGCAAGCGCTCTCCCAGCTGAGCTATAGCCCCGTCACGTAAAGCTTGTCGAGTTGACGGGCGGCA‑TCATATGAATCCCGCCCGAATGTGTCAACGGCAAATTGCaa                                           <  2:230/150‑1 (MQ=255)
                                                                                                            ggtggAGCTAAGCGGGATCGAACCGCTGACCTCTTGCATGCCATGCAAGCGCTCTCCCAGCTGAGCTATAGCCCCGTCACGTAAAGCTTGTCGAGTTGACGGGCGGCA‑TCATATGAATTCCGCCCGAATGTGTCAACGGCAAATTGCAAc                                          <  1:1468163/150‑1 (MQ=255)
                                                                                                             gtggAGCTAAGCGGGATCGAACCGCTGACCTCTTGCATGCCATGCAAGCGCTCTCCCAGCTGAGCTATAGCCCCGTCACGTAAAGCTTGTCGAGTTGACGGGCGGCA‑TCATATGAATTCCGCCCGAATGTGTCAACGGCAAATTGCAACg                                         <  2:711616/150‑1 (MQ=255)
                                                                                                             gtggAGCTAAGCGGGATCGAACCGCTGACCTCTTGCATGCCATGCAAGCGCTCTCCCAGCTGAGCTATAGCCCCGTCACGTAAAGCTTGTCGAGTTGACGGGCGGCA‑TCATATGAATTCCGCCCGAATGTGTCAACGGCAAATTGCAACg                                         <  2:1354298/150‑1 (MQ=255)
                                                                                                              tggAGCTAAGCGGGATCGAACCGCTGACCTCTTGCATGCCATGCAAGCGCTCTCCCAGCTGAGCTATAGCCCCGTCACGTAAAGCTTGTCGAGTTGACGGGCGGca                                                                                     <  1:1005301/106‑1 (MQ=255)
                                                                                                              tggAGCTAAGCGGGATCGAACCGCTGACCTCTTGCATGCCATGCAAGCGCTCTCCCAGCTGAGCTATAGCCCCGTCACGTAAAGCTTGTCGAGTTGACGGGCGGca                                                                                     >  2:1005301/1‑106 (MQ=255)
                                                                                                               ggAGCTAAGCGGGATCGAACCGCTGACCTCTTGCATGCCATGCAAGCGCTCTCCCAGCTGAGCTATAGCCCCGTCACGTAAAGCTTGTCGAGTTGACGGGCGGCA‑TCATATGAATTCCGCCCGAATGTGTCAACGGCAAATTGCAACgtg                                       >  1:31218/1‑150 (MQ=255)
                                                                                                                   cTAAGCGGGATCGAACCGCTGACCTCTTGCATGCCATGCAAGCGCTCTCCCAGCTGAGCTATAGCCCCGTCACGTAAAGCTTGTCGAGTTGACGGGCGGCA‑TCATATGAATTCCGCCCGAATGTGTCAACGGCAAATTGCAACGTGTAGt                                   >  2:337442/1‑150 (MQ=255)
                                                                                                                    tAAGCGGGATCGAACCGCTGACCTCTTGCATGCCATGCAAGCGCTCTCCCAGCTGAGCTATAGCCCCGTCACGTAAAGCTTGTCGAGTTGACGGGCGGCA‑TCATATGAATTCCGCCCGAATGTGTCAACGGCAAATTGCAACGTGTa                                     <  1:1332769/147‑1 (MQ=255)
                                                                                                                    tAAGCGGGATCGAACCGCTGACCTCTTGCATGCCATGCAAGCGCTCTCCCAGCTGAGCTATAGCCCCGTCACGTAAAGCTTGTCGAGTTGACGGGCGGCA‑TCATATGAATTCCGCCCGAATGTGTCAACGGCAAATTGCAACGTGTa                                     >  2:1332769/1‑147 (MQ=255)
                                                                                                                    tAAGCGGGATCGAACCGCTGACCTCTTGCATGCCATGCAAGCGCTCTCCCAGCTGAGCTATAGCCCCGTCACGTAAAGCTTGTCGAGTTGACGGGCGGCA‑TCATATGAATTCCGCCCGAATGTGTCAACGGCAAATTGCAACGTGTAGtt                                  >  1:906647/1‑150 (MQ=255)
                                                                                                                    tAAGCGGGATCGAACCGCTGACCTCTTGCATGCCATGCAAGCGCTCTCCCAGCTGAGCTATAGCCCCGTCACGTAAAGCTTGTCGAGTTGACGGGCGGCA‑TCATATGAATTCCGCCCGAATGTGTCAACGGCAAATTGCAACGTGTAGt                                   >  1:1339551/1‑149 (MQ=255)
                                                                                                                    tAAGCGGGATCGAACCGCTGACCTCTTGCATGCCATGCAAGCGCTCTCCCAGCTGAGCTATAGCCCCGTCACGTAAAGCTTGTCGAGTTGACGGGCGGCA‑TCATATGAATTCCGCCCGAATGTGTCAACGGCAAATTGCAACGTGTAGt                                   <  2:1339551/149‑1 (MQ=255)
                                                                                                                    tAAGCGGGATCGAACCGCTGACCTCTTGCATGCCATGCAAGCGCTCTCCCAGCTGAGCTATAGCCCCGTCACGTAAAGCTTGTCGAGTTGACGGCCGGCA‑TCATATGAATTCCGCACGAATGTGTCAACGGCAAATTGCAACGTGTAGtt                                  >  2:1664954/1‑150 (MQ=255)
                                                                                                                     aaGCGGGATCGAACCGCTGACCTCTTGCATGCCATGCAAGCGCTCTCCCAGCTGAGCTATAGCCCCGTCACGTAAAGCTTGTCGAGTTGACGGGCGGCA‑TCATATGAATTCCGCCCGAATGTGTCAACGGCAAATTGCAACGTGTAGttt                                 >  1:963473/1‑150 (MQ=255)
                                                                                                                      aGCGGGATCGAACCGCTGACCTCTTGCATGCCATGCAAGCGCTCTCCCAGCTGAGCTATAGCCCCGTCACGTAAAGCTTGTCGAGTTGACGGGCGGca‑t                                                                                   >  2:1389894/1‑99 (MQ=255)
                                                                                                                      aGCGGGATCGAACCGCTGACCTCTTGCATGCCATGCAAGCGCTCTCCCAGCTGAGCTATAGCCCCGTCACGTAAAGCTTGTCGAGTTGACGGGCGGca‑t                                                                                   <  1:1389894/99‑1 (MQ=255)
                                                                                                                       gCGGGATCGAACCGCTGACCTCTTGCATGCCATGCAAGCGCTCTCCCAGCTGAGCTATAGCCCCGTCACGTAAAGCTTGTCGAGTTGACGGGCGGCA‑TCATATGAATTCCGCCCGAATGTGTCAACGGCAAATTGCAACGTGTAGTTTCa                               >  1:1111545/1‑150 (MQ=255)
                                                                                                                        cGGGATCGAACCGCTGACCTCTTGCATGCCATGCAAGCGCTCTCCCAGCTGAGCTATAGCCCCGTCACGTAAAGCTTGTCGAGTTGACGGGCGGCA‑TCATATGAATTCCGCCCGAAtgtg                                                            <  2:1225746/120‑1 (MQ=255)
                                                                                                                        cGGGATCGAACCGCTGACCTCTTGCATGCCATGCAAGCGCTCTCCCAGCTGAGCTATAGCCCCGTCACGTAAAGCTTGTCGAGTTGACGGGCGGCA‑TCATATGAATTCCGCCCGAAtgtg                                                            >  1:1225746/1‑120 (MQ=255)
                                                                                                                          ggATCGAACCGCTGACCTCTTGCATGCCATGCAAGCGCTCTCCCAGCTGAGCTATAGCCCCGTCACGTAAAGCTTGTCGAGTTGACGGGCGGCA‑TCATATGAATTCCGCCCGAATGTGTCAACGGCAAATTGCAACGTGTAGTTTCAATc                            >  2:899688/1‑150 (MQ=255)
                                                                                                                          ggATCGAACCGCTGACCTCTTGCATGCCATGCAAGCGCTCTCCCAGCTGAGCTATAGCCCCGTCACGTAAAGCTTGTCGAGTTGACGGGCGGCA‑TCATATGAATTCCGCCCGAATGTGTCAACGGCAAATTGCAACGTGTAGTTTCAATc                            >  2:1441245/1‑150 (MQ=255)
                                                                                                                          ggATCGAACCGCTGACCTCTTGCATGCCATGCAAGCGCTCTCCCAGCTGAGCTATAGCCCCGTCACGTAAAGCTTGTCGAGTTGACGGGCGGCA‑TCATATGAATTCCGCCCGAATGTGTCAACGGCAAATTGCAACGTGTAGTTTCAATc                            <  1:186763/150‑1 (MQ=255)
                                                                                                                          ggATCGAACCGCTGACCTCTTGCATGCCATGCAAGCGCTCTCCCAGCTGAGCTATAGCCCCGTCAAGAAAAGCTTTTCCAGTTTACGGGCGGCA‑TAATATTAATTCCGCAAGAATGTGTCAACCGCAAAAACAAAAGAGTAGATTCAATc                            >  2:186763/1‑150 (MQ=18)
                                                                                                                           gATCGAACCGCTGACCTCTTGCATGCCATGCAAGCGCTCTCCCAGCTGAGCTATAGCCCCGTCACGTAAATCTTGTCGAGTTGACGGGCGGCA‑TCATATGAATTCCGCACCAATGTGTCAACGGCAAAATGCAAAGTGTAGTTTCAATCg                           >  2:1648130/1‑150 (MQ=37)
                                                                                                                           gATCGAACCGCTGACCTCTTGCATGCCATGCAAGCGCTCTCCCAGCTGAGCTATAGCCCCGTCACGTAAAGCTTGTCGAGTTGACGGGCGGCA‑TCATATGAATTCCGCCCGAATGTGTCAACGGCAAATTGCAACGTGTAGTTTCAATCg                           >  1:1312281/1‑150 (MQ=255)
                                                                                                                             tCGAACCGCTGACCTCTTGCATGCCATGCAAGCGCTCTCCCAGCTGAGCTATAGCCCCGTCACGTAAAGCTTGTCGAGTTGACGGGCGGCA‑TCATATGAATTCCGCCCGAATGTGTCAACGGCAAATTGCAACGTGTAGTTTCAATCGCt                         <  1:1155025/150‑1 (MQ=255)
                                                                                                                                  ccGCTGACCTCTTGCATGCCATGCAAGCGCTCTCCCAGCTGAGCTATAGCCCCGTCACGTAAAGCTTGTCGAGTTGACGGGCGGCA‑TCATATGAATTCCGCCCGAATGTGTCAACGGCAAATTGCAACGTGTAGTTTCAATCGCTGaaaa                    >  2:349060/1‑150 (MQ=255)
                                                                                                                                     cTGACCTCTTGCATGCCATGCAAGCGCTCTCCCAGCTGAGCTATAGCCCCGTCACGTAAAGCTTGTCGAGTTGACGGGCGGCA‑TCATATGAAATCCGCCCGAATGTGTCAACCGCAAATTGCAACGTGTAGTTTCAATCGCTGAAAAATc                 <  1:79174/150‑1 (MQ=255)
                                                                                                                                      tGACCTCTTGCATGCCATGCAAGCGCTCTCCCAGCTGAGCTATAGCCCCGTCACGTAAAGCTTGTCGAGTTGACGGGCGGCA‑TCATATGAATTCCGCCCGAATGTGTCAACGGCAAATTGCAACGTGTAGTTTCAATCGCTGAAAAATCa                <  1:1671077/150‑1 (MQ=255)
                                                                                                                                         gcgctTGCATGCCATGCAAGCGCTCTCCCAGCTGAGCTATAGCCCCGTCACGTAAAGCTTGTCGAGTTGACGGGCGGCA‑TCATATGAATTCCGCCCGAATGTGTCAACGGCAAATTGCAACGTGTAGTTTCAATCGCTGAAAAATTAGGc             >  1:1397011/4‑150 (MQ=255)
                                                                                                                                         ccTCTTGCATGCCATGCAAGCGCTCTCCCAGCTGAGCTATAGCCCCGTCACGTAAAGCTTGTCGAGTTGACGGGCGGCA‑TCATATGAATTCCGCGCGAATGTGTCAACGGCAAATTGCAACGTGTAGTTTCAATCGCTGAAAAATCAGGc             >  2:420488/1‑150 (MQ=255)
                                                                                                                                          cgctTGCATGGCATGCAAGCGCTCTCCCAGCTGAGCTATAGCCCCGTCACGTAAAGCTTGTCGAGTTGACGGGCGGCA‑TCATATGAATTCCGCCCGAATGTGTCAACGGCAAATTGCAACGTGTAGTTTCAATCGCTGAAAAATCAGGCa            >  1:457/3‑150 (MQ=255)
                                                                                                                                          cgctTGCATGGCATGCAAGCGCTCTCCCAGCTGAGCTATAGCCCCGTCACGTAAAGCTTGTCGAGTTGACGGGCGGCA‑TCATATGAATTCCGCCCGAATGTGTCAACGGCAAATTGCAACGTGTAGTTTCAATCGCTGAAAAATCAGGCa            >  1:415114/3‑150 (MQ=255)
                                                                                                                                           tctTGCATGCCATGCAAGCGCTCTCCCAGCTGAGCTATAGCCCCGTCACGTAAAGCTTGTCGAGTTGACGGGCGGCA‑TCATATGAATTCCGCCCGAATGTGTCAACGGCAAATTGCAACGTGTAGTTTCAATCGCTGAAAAATCAGGCaa           >  2:1158796/1‑150 (MQ=255)
                                                                                                                                             ttGCATGCCATGCAAGCGCTCTCCCAGCTGAGCTATAGCCCCGTCACGTAAAGCTTGTCGAGTTGACGGGCGGCA‑TCATATGAATTCCGCCCGAATGTGTCAACGGCAAATTGCAACGTGTAGTTTCAATCGCTGAAAAATCAGGCAAAt         >  2:323859/1‑150 (MQ=255)
                                                                                                                                              tGCATGCCATGCAAGCGCTCTCCCAGCTGAGCTATAGCCCCGTCACGTAAAGCTTGTCGAGTTGACGGGCGGCA‑TCATATGAATTCCGCCCGAATGTGTCAACGGCAAATTGCAACGTGTAGTTTCAATCGCTGAAAAATCAGGCAAATg        >  2:596157/1‑150 (MQ=255)
                                                                                                                                               gCATGCCATGCAAGCGCTCTCCCAGCTGAGCTATAGCCCCGTCACGTAAAGCTTGTCGAGTTGACGGGCGGCA‑TCATATGAATTCCGCCCGAATGTGTCAACGGCAAATTGCAACGTGTAGTTTCAATAGCTGAAAAATCAGGCAAATGa       >  2:1666002/1‑150 (MQ=255)
                                                                                                                                                catgccatgcAAGCGCTCTCCCAGCTGAGCTATAGCCCCGTCACGTAAAGCTTGTCGAGTTGACGGGCGGCA‑TCATATGAATTCCGCCCGAATGTGTCAACGGCAAATTGCAACGTGTAGTTTCAATCGCTGAAAAATCAGGCAAATGaa      <  2:391645/150‑1 (MQ=255)
                                                                                                                                                catgccatgcAAGCGCTCTCCCAGCTGAGCTATAGCCCCGTCACGTAAAGCTTGTCGAGTTGACGGGCGGCA‑TCATATGAATTCCGCCCGAATGTGTCAACGGCAAATTGCAACGTGTAGTTTCAATCGCTGAAAAATCAGGCAAATGaa      >  2:1043252/1‑150 (MQ=255)
                                                                                                                                                catgccatgcAAGCGCTCTCCCAGCTGAGCTATAGCCCCGTCACGTAAAGCTTGTCGAGTTGACGGGCGGCA‑TCATATGAATTCCGCCCGAATGTGTCAACGGCAAATTGCAACGTGTAGTTTCAATCGCTGAAAAATCAGGCAAATGaa      <  1:285413/150‑1 (MQ=255)
                                                                                                                                                 atgccatgcAAGCGCTCTCCCAGCTGAGCTATAGCCCCGTCACGTAAAGCTTGTCGAGTTGACGGGCGGCA‑TCATATGAATTCCGCCCGAATGTGTCAACGGCAAATTGCAACGTGTAGTTTCAATCGCTGAAAAATCAGGCAAATGAAc     <  2:979809/150‑1 (MQ=255)
                                                                                                                                                 atgccatgcAAGCGCTCTCCCAGCTGAGCTATAGCCCCGTCACGTAAAGCTTGTCGAGTTGACGGGCGGCA‑TCATATGAATTCCGCCCGAATGTGTCAACGGCAAATTGCAACGTGTAGTTTCAATCGCTGAAAAATCAGGCAAATGAAc     >  1:832493/1‑150 (MQ=255)
                                                                                                                                                 atgccatgcAAGCGCTCTCCCAGCTGAGCTATAGCCCCGTCACGTAAAGCTTGTCGAGTTGACGGGCGGCA‑TCATATGAATTCCGCCCGAATGTGTCAACGGCAAATTGCAACGTGTAGTTTCAATCGCTGAAAAATCAGGCAAATGAAc     >  1:4758/1‑150 (MQ=255)
                                                                                                                                                 atgccatgcAAGCGCTCTCCCAGCTGAGCTATAGCCCCGTCACGTAAAGCTTGTCGAGTTGACGGGCGGCA‑TCATATGAATTCCGCCCGAATGTGTCAACGGCAAATTGCAACGTGTAGTTTCAATCGCTGAAAAATCAGGCAAATGAAc     <  1:1648130/150‑1 (MQ=255)
                                                                                                                                                  tgccatgcAAGCGCTCTCCCAGCTGAGCTATAGCCCCGTCACGTAAAGCTTGTCGAGTTGACGGGCGGCA‑TCATATGAATTCCGCCCGAATGTGTCAACGGCAAATTGCAACGTGTAGTTTCAATCGCTGAAAAATCAGGCAAATGAACa    <  1:949878/150‑1 (MQ=255)
                                                                                                                                                    ccatgcAAGCGCTCTCCCAGCTGAGCTATAGCCCCGTCACGTAAAGCTTGTCGAGTTGACGGGCGGCA‑TCATATGAATTCCGCCCGAATGTGTCAACGGCAAATTGCAACGTGTAGTTTCAATCGCTGAAAAATCAGGCAAATGAACAAt  >  1:244682/1‑150 (MQ=255)
                                                                                                                                                    |                                                                                                                                                      
CTAAGCGGGATCGAACCGCTGACCTCTTGCATGCCATGCAAGCGCTCTCCCAGCTGAGCTATAGCCCCACGATGCGTTTACGTACCAAGTTTGCTGGGTGCAAAATTTGGTGGAGCTAAGCGGGATCGAACCGCTGACCTCTTGCATGCCATGCAAGCGCTCTCCCAGCTGAGCTATAGCCCCGTCACGTAAAGCTTGTCGAGTTGACGGGCGGCA‑TCATATGAATTCCGCCCGAATGTGTCAACGGCAAATTGCAACGTGTAGTTTCAATCGCTGAAAAATCAGGCAAATGAACAAT  >  NC_000913/2518049‑2518346

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: ‑