Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_000913 1,948,180 T→A 6.8% *151C (TGA→TGT nudB ← dihydroneopterin triphosphate pyrophosphatase

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0009131,948,1800TA 6.8% 263.8 / 8.1 74*151C (TGA→TGTnudBdihydroneopterin triphosphate pyrophosphatase
Reads supporting (aligned to +/- strand):  ref base T (33/36);  new base A (2/3);  total (35/39)
Fisher's exact test for biased strand distribution p-value = 1.00e+00
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 6.81e-01

CGCCCAGCTTAGCCGCGGTTACCAGCTCACGAATGATTTTAGTGAAGATTTTACCGCGCTTAGCATCCTGCGCAGCTTTACGATGTCTGGTGTTGGCCCATTTACTATGACCTGCCATAAAAATATCTCCAGATAGCCCTGCCTGTTCAGGCAGCGTTAATTACAAACTGTTCAATCGCCTGCCGGTTGCTCCAGGACTTAGTGAGCGCCGCCGCAGCAGGCGCATCAAGCCACTTGTAAGCCAGATGTTCAGTGAAAACGATCTGCCGCTCGTGCGGAAGCGCAAGACAGAA  >  NC_000913/1948034‑1948326
                                                                                                                                                  |                                                                                                                                                  
cGCCCAGCTTAGCCGCGGTTACCAGCTCACGAATGATTTTAGTGAAGATTTTACCGCGCTTAGCATCCTGCGCAGCTTTACGATGTCTGGTGTTGGCCCATTTACTATGACCTGCCATAAAAATATCTCCAGATAGCCCTGCCTGTTCa                                                                                                                                                  <  2:211123/149‑1 (MQ=255)
   ccAGCTTAGCCGCGGTTACCAGCTCACGAATGATTTTAGTGAAGATTTTACCGCGCTTAGCATCCTGCGCAGCTTTACGATGTCTGGTGTTGGCCCATTTACTATGACCTGCCATAAAAATATCTCCAGATAGCCCTGCCTGTTCa                                                                                                                                                  <  3:1333257/146‑1 (MQ=255)
   ccAGCTTAGCCGCGGTTACCAGCTCACGAATGATTTTAGTGAAGATTTTACCGCGCTTAGCATCCTGCGCAGCTTTACGATGTCTGGTGTTGGCCCATTTACTATGACCTGCCATAAAAATATCTCCAGATAGCCCTGCCTGAACAGGc                                                                                                                                               <  2:1200121/149‑1 (MQ=255)
   ccAGCTTAGCCGCGGTTACCAGCTCACGAATGATTTTAGCGAAGATTTTACCGCGCTTAGCATCCTGCGCAGCTTTACGATGTCTGGTGTTGGCCCATTTACTATGACCTGCCATAAAAATATCTCCAGATAGCCCTGCCTGTTCa                                                                                                                                                  >  4:1333257/1‑146 (MQ=255)
    cAGCTTAGCCGCGGTTACCAGCTCACGAATGATTTTAGTGAAGATTTTACCGCGCTTAGCATCCTGCGCAGCTTTACGATGTCTGGTGTTGGCCCATTTACTATGACCTGCCATAAAAATATCTCCAGATAGCCCTGCCTGTTc                                                                                                                                                   <  1:1206242/144‑1 (MQ=255)
    cAGCTTAGCCGCGGTTACCAGCTCACGAATGATTTTAGTGAAGATTTTACCGCGCTTAGCATCCTGCGCAGCTTTACGATGTCTGGTGTTGGCCCATTTACTATGACCTGCCATAAAAATATCTCCAGATAGCCCTGCCTGTTc                                                                                                                                                   >  2:1206242/1‑144 (MQ=255)
    cAGCTTAGCCGCGGTTACCAGCTCACGAATGATTTTAGTGAAGATTTTACCGCGCTTAGCATCCTGCGCAGCTTTACGATGTCTGGTGTTGGCCCATTTACTATGACCTGCCATAAAAATATCTCCAGATAGCCCTGCCTGTTCAGGCa                                                                                                                                              <  3:399816/149‑1 (MQ=255)
           gCCGCGGTTACCAGCTCACGAATGATTTTAGTGAAGATTTTACCGCGCTTAGCATCCTGCGCAGCTTTACGATGTCTGGTGTTGGCCCATTTACTATGACCTGCCATAAAAATATCTCCAGATAGCCCTGCCTGTTCAGGCAGCGTTaa                                                                                                                                       <  1:1158898/149‑1 (MQ=255)
           gCCGCGGTTACCAGCTCACGAATGATTTTAGTGAAGATTTTACCGCGCCTAGCATCCTGCGCAGCTTTACGATGTCTGGTGTTGGCCCATTTACTATGACCTGCCATAAAAATATCTCCAGATAGCCCTGCCTGTTCAGGCAGCGTTaa                                                                                                                                       <  3:1209915/149‑1 (MQ=255)
            ccGCGGTTACCAGCTCACGAATGATTTTAGTGAAGATTTTACCGCGCTTAGCATCCTGCGCAGCTTTACGATGTCTGGTGTTGGCCCATTTACTATGACCTGCCATAAAAATATCTCCAGATAGCCCTGCCTGTTCAg                                                                                                                                                 <  1:1417941/138‑1 (MQ=255)
            ccGCGGTTACCAGCTCACGAATGATTTTAGTGAAGATTTTACCGCGCTTAGCATCCTGCGCAGCTTTACGATGTCTGGTGTTGGCCCATTTACTATGACCTGCCATAAAAATATCTCCAGATAGCCCTGCCTGTTCAg                                                                                                                                                 >  2:1417941/1‑138 (MQ=255)
                   tACCAGCTCACGAATGATTTTAGTGAAGATCTTACCGCGCTTAGCATCCTGCGCAGCTTTACGATGTCTGGTGTTGGCCCATTTACTATGACCTGCCATAAAAATATCTCCAGATAGCCCTGCCTGTTCAGGCAGCGTTAATTACAAAc                                                                                                                               >  4:1315411/1‑149 (MQ=255)
                    aCCAGCTCACGAATGATTTTAGTGAAGATTTTACCGCGCTTAGCATCCTGCGCAGCTTTACGATGTCTGGTGTTGGCCCATTTACTATGACCTGCCATAAAAATATCTCCAGATAGCCCTGCCTGTTCAGGCAGCGTTAAt                                                                                                                                      <  3:1305640/141‑1 (MQ=255)
                    aCCAGCTCACGAATGATTTTAGTGAAGATTTTACCGCGCTTAGCATCCTGCGCAGCTTTACGATGTCTGGTGTTGGCCCATTTACTATGACCTGCCATAAAAATATCTCCAGATAGCCCTGCCTGTTCAGGCAGCGTTAAt                                                                                                                                      >  4:1305640/1‑141 (MQ=255)
                       aGCTCACGAATGATTTTAGTGAAGATTTTACCGCGCTTAGCATCCTGCGCAGCTTTACGATGTCTGGTGTTGGCCCATTTACTATGACCTGCCATAAAAATATCTCCAGATAGCCCTGCCTGTTCAGGCAGCGTTAAt                                                                                                                                      <  1:1449029/138‑1 (MQ=255)
                       aGCTCACGAATGATTTTAGTGAAGATTTTACCGCGCTTAGCATCCTGCGCAGCTTTACGATGTCTGGTGTTGGCCCATTTACTATGACCTGCCATAAAAATATCTCCAGATAGCCCTGCCTGTTCAGGCAGCGTTAAt                                                                                                                                      >  2:1449029/1‑138 (MQ=255)
                            aCGAATGATTTTAGTGAAGATTTTACCGCGCTTAGCATCCTGCGCAGCTTTACGATGTCTGGTGTTGGCCCATTTACTATGACCTGCCATAAAAATATCTCCAGATAGCCCTGCCTGTTCAGGCAGCGTTAATTACAAACTGTTCAATc                                                                                                                      <  2:134038/149‑1 (MQ=255)
                                aTGATTTTAGTGAAGATCTTACCGCGCTTAGCATCCTGCGCAGCTTTACGATGTCTGGTGTTGGCCCATTTACTATGACCTGCCATAAAAATATCTCCAGATAGCCCTGCCTGTTCAGGCAGCGTTAATTACAAACTGTTCAATCGCCt                                                                                                                  <  3:1315411/149‑1 (MQ=255)
                                            aaGATTTTACCGCGCTTAGCATCCTGCGCAGCTTTACGATGTCTGGTGTTGGCCCATTTACTATGACCTGCCATAAAAATATCTCCAGATAGCCCTGCCTGTTCAGGCAGCGTTAATTACAAACTGTTCAATCGCCTGCCGGTTGCTcc                                                                                                      >  4:267699/1‑149 (MQ=255)
                                              gATTTTACCGCGCTTAGCATCCTGCGCAGCTTTACGATGTCTGGTGTTGGCCCATTTACTATGACCTGCCATAAAAATATCTCCAGATAGCCCTGCCTGAACAGGCAggg                                                                                                                                           <  1:1433722/110‑3 (MQ=255)
                                              gATTTTACCGCGCTTAGCATCCTGCGCAGCTTTACGATGTCTGGTGTTGGCCCATTTACTATGACCTGCCATAAAAATATCTCCAGATAGCCCTGCCTGAACAGGCAggg                                                                                                                                           >  2:1433722/1‑108 (MQ=255)
                                                ttttACCGCGCTTAGCATCCTGCGCAGCTTTACGATGTCTGGTGTTGGCCCATTTACTATGACCTGCCATAAAAATATCTCCAGATAGCCCTGCCTGTTCAGGCAGCGTTAATTACAAACTGTTCAATCGCCTGCCGGTTGCTCCAg                                                                                                    <  2:863751/147‑1 (MQ=255)
                                                ttttACCGCGCTTAGCATCCTGCGCAGCTTTACGATGTCTGGTGTTGGCCCATTTACTATGACCTGCCATAAAAATATCTCCAGATAGCCCTGCCTGTTCAGGCAGCGTTAATTACAAACTGTTCAATCGCCTGCCGGTTGCTCCAg                                                                                                    >  1:863751/1‑147 (MQ=255)
                                                    aCCGCGCTTAGCATCCTGCGCAGCTTTACGATGTCTGGTGTTGGCCCATTTACTATGACCTGCCATAAAAATATCTCCAGATAGCCCTGCCTGTTCAGGCAGCGTTAATTACAAACTGTTCAATc                                                                                                                      <  3:383948/125‑1 (MQ=255)
                                                    aCCGCGCTTAGCATCCTGCGCAGCTTTACGATGTCTGGTGTTGGCCCATTTACTATGACCTGCCATAAAAATATCTCCAGATAGCCCTGCCTGTTCAGGCAGCGTTAATTACAAACTGTTCAATc                                                                                                                      >  4:383948/1‑125 (MQ=255)
                                                          cTTAGCATCCTGCGCAGCTTTACGATGTCTGGTGTTGGCCCATTTACTATGACCTGCCATAAAAATATCTCCAGATAGCCCTGCCTGTTCAGGCAGCGTTAATTACAAACTGTTCAATCGCCTGCCGGTTGCTCCAg                                                                                                    <  3:218660/137‑1 (MQ=255)
                                                          cTTAGCATCCTGCGCAGCTTTACGATGTCTGGTGTTGGCCCATTTACTATGACCTGCCATAAAAATATCTCCAGATAGCCCTGCCTGTTCAGGCAGCGTTAATTACAAACTGTTCAATCGCCTGCCGGTTGCTCCAg                                                                                                    >  4:218660/1‑137 (MQ=255)
                                                              gCATCCTGCGCAGCTTTACGATGTCTGGTGTTGGCCCATTTACTATGACCTGCCATAAAAATATCTCCAGATAGCCCTGCCTGAACAGGCAggg                                                                                                                                           <  1:40742/94‑3 (MQ=255)
                                                              gCATCCTGCGCAGCTTTACGATGTCTGGTGTTGGCCCATTTACTATGACCTGCCATAAAAATATCTCCAGATAGCCCTGCCTGAACAGGCAggg                                                                                                                                           >  2:40742/1‑92 (MQ=255)
                                                               cATCCTGCGCAGCTTTACGATGTCTGGTGTTGGCCCATTTACTATGACCTGCCATAAAAATATCTCCAGATAGCCCTGCCTGTTCAGGCAGCGTTAATTACAAACTGTTCAATCGCCTGCCGGTTGCTCCAGGACTTa                                                                                              <  2:543461/138‑1 (MQ=255)
                                                               cATCCTGCGCAGCTTTACGATGTCTGGTGTTGGCCCATTTACTATGACCTGCCATAAAAATATCTCCAGATAGCCCTGCCTGTTCAGGCAGCGTTAATTACAAACTGTTCAATCGCCTGCCGGTTGCTCCAGGACTTa                                                                                              >  1:543461/1‑138 (MQ=255)
                                                                  ccTGCGCAGCTTTACGATGTCTGGTGTTGGCCCATTTACTATGACCTGCCATAAAAATATCTCCAGATAGCCCTGCCTGTTCAGGCAGCGTTAATTACAAACTGTTCAATCGCCTGCCGGTTGCTCCAGGACTTAGTGAGcgccgccgc                                                                                <  3:267699/149‑1 (MQ=255)
                                                                   cTGCGCAGCTTTACGATGTCTGGTGTTGGCCCATTTACTATGACCTGCCATAAAAATATCTCCAGATAGCCCTGCCTGTTCAGGCAGCGTTAATTACAAACTGTTCAATCGCCTGCCGGTTGCTCCAGGACTTAGTg                                                                                           >  3:293631/1‑137 (MQ=255)
                                                                   cTGCGCAGCTTTACGATGTCTGGTGTTGGCCCATTTACTATGACCTGCCATAAAAATATCTCCAGATAGCCCTGCCTGTTCAGGCAGCGTTAATTACAAACTGTTCAATCGCCTGCCGGTTGCTCCAGGACTTAGTg                                                                                           <  4:293631/137‑1 (MQ=255)
                                                                    tGCGCAGCTTTACGATGTCTGGTGTTGGCCCATTTACTATGACCTGCCATAAAAATATCTCCAGATAGCCCTGCCTGTTCAGGCAGc                                                                                                                                            >  3:188757/1‑87 (MQ=255)
                                                                    tGCGCAGCTTTACGATGTCTGGTGTTGGCCCATTTACTATGACCTGCCATAAAAATATCTCCAGATAGCCCTGCCTGTTCAGGCAGc                                                                                                                                            <  4:188757/87‑1 (MQ=255)
                                                                         aGCTTTACGATGTCTGGTGTTGGCCCATTTACTATGACCTGCCATAAAAATATCTCCAGATAGCCCTGCCTGTTCAGGCAGCGTTAATTACAAACTGTTCAATCGCCTGCCGGTTGCTCCAGGACTTAGTGAGCGCCGCCGCAGCAGgc                                                                         >  2:505912/1‑149 (MQ=255)
                                                                          gCTTTACGATGTCTGGTGTTGGCCCATTTACTATGACCTGCCATAAAAATATCTCCAGATAGCCCTGCCTGTTCAGGCAGCGTTAATTACAAACTGTTCAATCGCCTGCCGGTTGCTCCAg                                                                                                    >  3:779051/1‑121 (MQ=255)
                                                                          gCTTTACGATGTCTGGTGTTGGCCCATTTACTATGACCTGCCATAAAAATATCTCCAGATAGCCCTGCCTGTTCAGGCAGCGTTAATTACAAACTGTTCAATCGCCTGCCGGTTGCTCCAg                                                                                                    <  4:779051/121‑1 (MQ=255)
                                                                            tttACGATGTCTGGTGTTGGCCCATTTACTATGACCTGCCATAAAAATATCTCCAGATAGCCCTGCCTGTTCAGGCAGCGTTAATTACAAACTGTTCAATCGCCTGCCGGTTGCTCCAGGACTTAGTGAGCGCCGCCgcag                                                                              >  2:601055/1‑141 (MQ=255)
                                                                            tttACGATGTCTGGTGTTGGCCCATTTACTATGACCTGCCATAAAAATATCTCCAGATAGCCCTGCCTGTTCAGGCAGCGTTAATTACAAACTGTTCAATCGCCTGCCGGTTGCTCCAGGACTTAGTGAGCGCCGCCgcag                                                                              <  1:601055/141‑1 (MQ=255)
                                                                                               gCCCATTTACTATGACCTGCCATAAAAATATCTCCAGATAGCCCTGCCTGTTCAGGCAGCGTTAATTACAAACTGTTCAATCGCCTGCCGGTTGCTCCAGGACTTAGTGAGCGCCGCCGCAGCAGGCGCATCaa                                                                  <  2:591700/134‑1 (MQ=255)
                                                                                               gCCCATTTACTATGACCTGCCATAAAAATATCTCCAGATAGCCCTGCCTGTTCAGGCAGCGTTAATTACAAACTGTTCAATCGCCTGCCGGTTGCTCCAGGACTTAGTGAGCGCCGCCGCAGCAGGCGCATCaa                                                                  >  1:591700/1‑134 (MQ=255)
                                                                                                   aTTTACTATGACCTGCCATAAAAATATCTCCAGATAGCCCTGCCTGTTCAGGCAGCGTTAATTACAAACTGTTCAATCGCCTGCCGGTTGCTCCAGGACTTAGTGAGCGCCGCCGCAGCAGGCGCATc                                                                    <  1:1191043/128‑1 (MQ=255)
                                                                                                   aTTTACTATGACCTGCCATAAAAATATCTCCAGATAGCCCTGCCTGTTCAGGCAGCGTTAATTACAAACTGTTCAATCGCCTGCCGGTTGCTCCAGGACTTAGTGAGCGCCGCCGCAGCAGGCGCATc                                                                    >  2:1191043/1‑128 (MQ=255)
                                                                                                         tATGACCTGCCATAAAAATATCTCCAGATAGCCCTGCCTGTTCAGGCAGCGTTAATTACAAACTGTTCAATCGCCTGCCGGTTGCTCCAGGACTTAGTGAGCGCCGCCGCAGCAGGCGCATCAAGCCACTTGTAAGcc                                                    <  1:1301598/138‑1 (MQ=255)
                                                                                                         tATGACCTGCCATAAAAATATCTCCAGATAGCCCTGCCTGTTCAGGCAGCGTTAATTACAAACTGTTCAATCGCCTGCCGGTTGCTCCAGGACTTAGTGAGCGCCGCCGCAGCAGGCGCATCAAGCCACTTGTAAGcc                                                    >  2:1301598/1‑138 (MQ=255)
                                                                                                            gACCTGCCATAAAAATATCTCCAGATAGCCCTGCCTGTTCAGGCAGCGTTAATTACAAACTGTTCAATCGCCTGCCGGTTGCTCCAGGACTTAGTGAGCGCCGCCGcagcag                                                                           >  4:261561/1‑112 (MQ=255)
                                                                                                            gACCTGCCATAAAAATATCTCCAGATAGCCCTGCCTGTTCAGGCAGCGTTAATTACAAACTGTTCAATCGCCTGCCGGTTGCTCCAGGACTTAGTGAGCGCCGCCGcagcag                                                                           <  3:261561/112‑1 (MQ=255)
                                                                                                                tGCCATAAAAATATCTCCAGATAGCCCTGCCTGTTCAGGCAGCGTTAATTACAAACTGTTCAATCGCCTGCCGGTTGCTCCAGGACTTa                                                                                              >  1:908250/1‑89 (MQ=255)
                                                                                                                tGCCATAAAAATATCTCCAGATAGCCCTGCCTGTTCAGGCAGCGTTAATTACAAACTGTTCAATCGCCTGCCGGTTGCTCCAGGACTTa                                                                                              >  2:1368054/1‑89 (MQ=255)
                                                                                                                tGCCATAAAAATATCTCCAGATAGCCCTGCCTGTTCAGGCAGCGTTAATTACAAACTGTTCAATCGCCTGCCGGTTGCTCCAGGACTTa                                                                                              <  1:1368054/89‑1 (MQ=255)
                                                                                                                tGCCATAAAAATATCTCCAGATAGCCCTGCCTGTTCAGGCAGCGTTAATTACAAACTGTTCAATCGCCTGCCGGTTGCTCCAGGACTTa                                                                                              <  2:908250/89‑1 (MQ=255)
                                                                                                                  ccATAAAAATATCTCCAGATAGCCCTGCCTGTTCAGGCAGCGTTAATTACAAACTGTTCAATCGCCTGCCGGTTGCTCCAGGACTTAGTGAGCGCCGCCGCAGCAGGCGCATCAAGCCACTTGTAAGCCAGATGTTCAGTGAAAACGAt                                <  1:505912/149‑1 (MQ=255)
                                                                                                                  ccATAAAAATATCTCCAGATAGCCCTGCCTGTTCAGGCAGCGTTAATTACAAACTGTTCAATCGCCTGCCGGTTGCTCCAGGACTTAGTGAGCGCCGCCGCAGCAGGCGCATCAAGCCACTTGTAAGCCAGATGTTCAGTGAAAACGAt                                >  4:147649/1‑149 (MQ=255)
                                                                                                                     tAAAAATATCTCCCGATAGCCCTGCCTGTTCAGGCAGCGTTAATTACAAACTGTTCAATCGCCTGCCGGTTGCTCCAGGACTTAGTGAGCGCCGCCGCAGCAGGCGCATCAAGCCACTTGt                                                         <  1:1027937/121‑1 (MQ=255)
                                                                                                                     tAAAAATATCTCCAGATAGCCCTGCCTGTTCAGGCAGCGTTAATTACAAACTGTTCAATCGCCTGCCGGTTGCTCCAGGACTTAGTGAGCGCCGCCGCAGCAGGCGCATCAAGCCACTTGt                                                         >  2:1027937/1‑121 (MQ=255)
                                                                                                                         aaTATCTCCAGATAGCCCTGCCTGTTCAGGCAGCGTTAATTACAAACTGTTCAATCGCCTGCCGGTTGCTCCAGGACTTAGTGAGCGCCGCCGCAGCAGGCGCATCAAGCCACTTGTAAGCCAGATGTTCAGTGAAAACGATCTGCCGc                         >  1:883418/1‑149 (MQ=255)
                                                                                                                               tcCAGATAGCCCTGCCTGTTCAGGCAGCGTTAATTACAAACTGTTCAATCGCCTGCCGGTTGCTCCAGGACTTAGTGAGCGCCGCCGCAGCAGGCGCATCAAGCCAc                                                             <  1:671805/107‑1 (MQ=255)
                                                                                                                               tcCAGATAGCCCTGCCTGTTCAGGCAGCGTTAATTACAAACTGTTCAATCGCCTGCCGGTTGCTCCAGGACTTAGTGAGCGCCGCCGCAGCAGGCGCATCAAGCCAc                                                             >  2:671805/1‑107 (MQ=255)
                                                                                                                               tcCAGATAGCCCTGCCTGTTCAGGCAGCGTTAATTACAAACTGTTCAATCGCCTGCCGGTTGCTCCAGGACTTAGTGAGCGCCGCCGCAGCAGGCGCATCAAGCCACTTGTAAGCCAGATGTTCAGTGAAAACGATCTGCCGCTCGTGc                   >  3:64612/1‑149 (MQ=255)
                                                                                                                               tcCAGATAGCCCTGCCTGTTCAGGCAGCGTTAATTACAAACTGTTCAATCGCCTGCCGGTTGCTCCAGGACTTAGTGAGCGCCGCCGCAGCAGGCGCATCAAGCCACTTGTAAGCCAGATGTTCAGTGAAAACGATCTGCCGCTCGTGc                   <  2:883418/149‑1 (MQ=255)
                                                                                                                                        cccTGCCTGTTCAGGCAGCGTTAATTACAAACTGTTCAATCGCCTGCCGGTTGCTCCAGGACTTAGTGAGCGCCGCCGCAGCAGGCGCATCAAGCCACTTGTAAGCCAGATGTTCAGTGAAAACGATCTGCCGCTCGTGCGGAAGCGCa          >  2:1404756/1‑149 (MQ=255)
                                                                                                                                        cccTGCCTGTTCAGGCAGCGTTAATTACAAACTGTTCAATCGCCTGCCGGTTGCTCCAGGACTTAGTGAGCGCCGCCGCAGCAGGCGCATCAAGCCACTTGTAAGCCAGATGTTCAGTGAAAACGATCTGCCGCTCGTGCGGAAGCGCa          >  1:1236921/1‑149 (MQ=255)
                                                                                                                                         cctgcctgTTCAGGCAGCGTTAATTACAAACTGTTCAATCGCCTGCCGGTTGCTCCAGGACTTAGTGAGCGCCGCCGCAGCAGGCGCATCAAGCCACTTGTAAGCCAGATGTTCAGTGAAAACGATCTGCCGCTCGTGCGGAAGCGCaa         <  2:1236921/149‑1 (MQ=255)
                                                                                                                                          ctgcctgTTCAGGCAGCGTTAATTACAAACTGTTCAATCGCCTGCCGGTTGCTCCAGGACtt                                                                                               <  3:590398/62‑1 (MQ=255)
                                                                                                                                          ctgcctgTTCAGGCAGCGTTAATTACAAACTGTTCAATCGCCTGCCGGTTGCTCCAGGACtt                                                                                               >  4:590398/1‑62 (MQ=255)
                                                                                                                                          ctgcctgTTCAGGCAGCGTTAATTACAAACTGTTCAATCGCCTGCCGGTTGCTCCAGGACTTAGTGAGCGCCGCCGCAGCAGGCGCATCAAGCCACTTGTAAGCCAGAt                                                >  4:1070002/1‑109 (MQ=255)
                                                                                                                                          ctgcctgTTCAGGCAGCGTTAATTACAAACTGTTCAATCGCCTGCCGGTTGCTCCAGGACTTAGTGAGCGCCGCCGCAGCAGGCGCATCAAGCCACTTGTAAGCCAGAt                                                <  3:1070002/109‑1 (MQ=255)
                                                                                                                                          ctgcctgTTCAGGCAGCGTTAATTACAAACTGTTCAATCGCCTGCCGGTTGCTCCAGGACTTAGTGAGCGCCGCCGCAGCAGGCGCATCAAGCCACTTGTAAGCCAGATGtt                                             <  1:741941/112‑1 (MQ=255)
                                                                                                                                          ctgcctgTTCAGGCAGCGTTAATTACAAACTGTTCAATCGCCTGCCGGTTGCTCCAGGACTTAGTGAGCGCCGCCGCAGCAGGCGCATCAAGCCACTTGTAAGCCAGATGtt                                             >  2:741941/1‑112 (MQ=255)
                                                                                                                                          ctgcctgTTCAGGCAGCGTTAATTACAAACTGTTCAATCGCCTGCCGGTTGCTCCAGGACTTAGTGAGCGCCGCCGCAGCAGGCGCATCAAGCCACTTGTAAGCCAGATGTTCAGTGAAAACGATCTGCCGCTCGTGCGGa                <  4:1375754/141‑1 (MQ=255)
                                                                                                                                          ctgcctgTTCAGGCAGCGTTAATTACAAACTGTTCAATCGCCTGCCGGTTGCTCCAGGACTTAGTGAGCGCCGCCGCAGCAGGCGCATCAAGCCACTTGTAAGCCAGATGTTCAGTGAAAACGATCTGCCGCTCGTGCGGa                >  3:1375754/1‑141 (MQ=255)
                                                                                                                                                gTTCAGGCAGCGTTAATTACAAACTGTTCAATCGCCTGCCGGTTGCTCCAGGACTTAGTGAGCGCCGCCGCAGCAGGCGCATCAAGCCACTTGTAAGCCAGATGTTCAGTGAAAACGATCTGCCGCTCGTGCGGAAGCGCAAGACAGaa  <  1:1404756/149‑1 (MQ=255)
                                                                                                                                                 ttCAGGCAGCGTTAATTACAAACTGTTCAATCGCCTGCCGGTTGCTCCAGGACTTAGTGAGCGCCGCCGCAGCAGGCGCATCAAGCCACTTGTAAGCCAGATGTTCAGTGAAAAc                                   >  2:587263/1‑115 (MQ=255)
                                                                                                                                                 ttCAGGCAGCGTTAATTACAAACTGTTCAATCGCCTGCCGGTTGCTCCAGGACTTAGTGAGCGCCGCCGCAGCAGGCGCATCAAGCCACTTGTAAGCCAGATGTTCAGTGAAAAc                                   <  1:587263/115‑1 (MQ=255)
                                                                                                                                                  |                                                                                                                                                  
CGCCCAGCTTAGCCGCGGTTACCAGCTCACGAATGATTTTAGTGAAGATTTTACCGCGCTTAGCATCCTGCGCAGCTTTACGATGTCTGGTGTTGGCCCATTTACTATGACCTGCCATAAAAATATCTCCAGATAGCCCTGCCTGTTCAGGCAGCGTTAATTACAAACTGTTCAATCGCCTGCCGGTTGCTCCAGGACTTAGTGAGCGCCGCCGCAGCAGGCGCATCAAGCCACTTGTAAGCCAGATGTTCAGTGAAAACGATCTGCCGCTCGTGCGGAAGCGCAAGACAGAA  >  NC_000913/1948034‑1948326

Alignment Legend
Aligned base mismatch/match (shaded by quality score): ATCG/ATCG < 3 ≤ ATCG/ATCG < 27 ≤ ATCG/ATCG < 34 ≤ ATCG/ATCG < 39 ≤ ATCG/ATCG < 40 ≤ ATCG/ATCG
Unaligned base: atcg    Masked matching base: atcg    Alignment gap:     Deleted base: