Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_000913 3,431,160 A→C 6.1% F86V (TTT→GTT)  tsaC ← threonylcarbamoyl‑AMP synthase

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0009133,431,1600AC 6.1% 312.9 / 10.2 98F86V (TTT→GTT) tsaCthreonylcarbamoyl‑AMP synthase
Reads supporting (aligned to +/- strand):  ref base A (51/41);  new base C (2/4);  total (53/45)
Fisher's exact test for biased strand distribution p-value = 4.09e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 9.95e-01

GGTTTACCATAAGCCTGGCACAAAGCAACCACCAACGGATGGTCGGTGACTCGTACAGCAAGCGAATCAAAGCGGCCCGTCAACCAGCGCGGTGTTGTCGCAGGCGCGGGAAAGACAAAGGTGACAGGACCTGGCCAGCGGGAAAAAATGGTTTCACGCTGCACGTCAGTCAACATGGTGTCATCAATATAGGGTTTAAGCTGCTCGTAATTTGCTGCGATTAAAATCAGCCCCTTATCAACCGGACGCTGTTTTAACTCCAACAGTCGCATCACTGCTGTTTCGCTATCAGGA  >  NC_000913/3431015‑3431308
                                                                                                                                                 |                                                                                                                                                    
ggTTTACCATAAGCCTGGCACAAAGCAACCACCAACGGATGGTCGGTGACTCGTACAGCAAGCGAATCAAAGCGGCCCGTCAACCAGCGCGGTGTTGTCGCAGGCGCGGGAAAGACAAAGGTGACAGGACCTGGCCAGCGGGAAAAAATg                                                                                                                                                  <  1:707529/150‑1 (MQ=255)
 gTTTACCATAAGCCTGGCACAAAGCAACCACCAACGGATGGTCGGTGACTCGTACAGCAAGCGAATCAAAGCGGCCCGTCAACCAGCGCGGTGTTGTCGCAGGCGCGGGAAAGACAAAGGTGACAGGACCTGGCCAGCGGGAAAAAATgg                                                                                                                                                 >  2:92142/1‑150 (MQ=255)
 gTTTACCATAAGCCTGGCACAAAGCAACCACCAACGGATGGTCGGTGACTCGTACAGCAAGCGAATCAAAGCGGCCCGTCAACCAGCGCGGTGTTGTCGCAGGCGCGGGAAAGACAAAGGTGACAGGACCTGGCCAGCGGGAAAAAATgg                                                                                                                                                 >  2:97748/1‑150 (MQ=255)
    tACCATAAGCCTGGCACAAAGCAACCACCAACGGATGGTCGGTGACTCGTACAGCAAGCGAATCAAAGCGGCCCGTCAACCAGCGCGGTGTTGTCGCAGGCGCGGGAAAGACAAAGGTGACAGGACCTGGCCAGCGGGAAAAAATGGttt                                                                                                                                              <  1:6144/150‑1 (MQ=255)
    tACCATAAGCCTGGCACAAAGCAACCACCAACGGATGGTCGGTGACTCGTACAGCAAGCGAATCAAAGCGGCCCGTCAACCAGCGCGGTGTTGTCGCAGGCGCGGGAAAGACAAAGGTGACAGGACCTGGCCAGCGGGAAAAAATGGttt                                                                                                                                              >  1:204652/1‑150 (MQ=255)
    tACCATAAGCCTGGCACAAAGCAACCACCAACGGATGGTCGGTGACTCGTACAGCAAGCGAATCAAAGCGGCCCGTCAACCAGCGCGGTGTTGTCGCAGGCGCGGGAAAGACAAAGGTGACAGGACCTGGCCAGCGGGAAAAAATGGttt                                                                                                                                              >  2:613589/1‑150 (MQ=255)
     aCCATAAGCCTGGCACAAAGCAACCACCAACGGATGGTCGGTGACTCGTACAGCAAGCGAATCAAAGCGGCCCGTCAACCAGCGCGGTGTTGTCGCAGGCGCGGGAAAGACAAAGGTGACAGGACCTGGCCAGCGGGAAAAAATGGTTTc                                                                                                                                             >  1:908710/1‑150 (MQ=255)
     aCCATAAGCCTGGCACAAAGCAACCACCAACGGATGGTCGGTGACTCGTACAGCAAGCGAATCAAAGCGGCCCGTCAACCAGCGCGGTGTTGTCGCAGGCGCGGGAAAGACAAAGGTGACAGGACCTGGCCAGCGGGAAAAAATGGTTTc                                                                                                                                             >  1:87527/1‑150 (MQ=255)
       cATAAGCCTGGCACAAAGCAACCACCAACGGATGGTCGGTGACTCGTACAGCAAGCGAATCAAAGCGGCCCGTCAACCAGCGCGGTGTTGTCGCAGGCGCGGGAAAGACAAAGGTGACAGGACCTGGCCAGCGGGAAAAAATGGTTTCAc                                                                                                                                           >  1:1268264/1‑150 (MQ=255)
       cATAAGCCTGGCACAAAGCAACCACCAACGGATGGTCGGTGACTCGTACAGCAAGCGAATCAAAGCGGCCCGTCAACCAGCGCGGTGTTGTCGCAGGCGCGGGAAAGACAAAGGTGACAGGACCTGGCCAGCGGGAAAAAATGGTTTCAc                                                                                                                                           <  2:204652/150‑1 (MQ=255)
          aaGCCTGGCACAAAGCAACCACCAACGGATGGTCGGTGACTCGTACAGCAAGCGAATCAAAGCGGCCCGTCAACCAGCGCGGTGTTGTCGCAGGCGCGGGAAAGACAAAGGTGACAGGACCTGGCCAGCGGGAAAAAATGGTTTCACGCt                                                                                                                                        >  2:1496166/1‑150 (MQ=255)
           aGCCTGGCACAAAGCAACCACCAACGGATGGTCGGTGACTCGTACAGCAAGCGAATCAAAGCGGCCCGTCAACCAGCGCGGTGTTGTCGCAGGCGCGGGAAAGACAAAGGTGACAGGACCTGGCCAGCGTGAAACCATTTTTTCCCGCTg                                                                                                                                       <  2:882466/150‑1 (MQ=255)
           aGCCTGGCACAAAGCAACCACCAACGGATGGTCGGTGACTCGTACAGCAAGCGAATCAAAGCGGCCCGTCAACCAGCGCGGTGTTGTCGCAGGCGCGGGAAAGACAAAGGTGACAGGACCTGGCCAGCGGGAAAAAATGGTTTCACGCTg                                                                                                                                       >  2:1132567/1‑150 (MQ=255)
           aGCCTGGCACAAAGCAACCACCAACGGATGGTCGGTGACTCGTACAGCAAGCGAATCAAAGCGGCCCGTCAACCAGCGCGGTGTTGTCGCAGGCGCGGGAAAGACAAAGGTGACAGGACCTGGCCAGCGGGAAAAAATGGTTTCACGCTg                                                                                                                                       <  2:364294/150‑1 (MQ=255)
           aGCCTGGCACAAAACAACCACCAACGGATGGTCGGTGACTCGTACAGCAAGCGAATCAAAGCGGCCCGTCAACCAGCGCGGTGTTGTCGCAGGCGCGGGAAAGACAAAGGTGACAGGACCTGGCCAGCGGGAAAAAATGGTTTCACGCTg                                                                                                                                       <  1:1393573/150‑1 (MQ=255)
            gCCTGGCACAAAGCAACCACCAACGGATGGTCGGTGACTCGTACAGCAAGCTAATCAAAGCGGCCCGTCAACCAGCGCTTTGTTTTCGCAGTCGCGTGAAATACAAAGGTGACAGGACCTGTCCAGCGTGAAACCATTTTTTCCCGCTgg                                                                                                                                      <  2:1063816/150‑2 (MQ=255)
             ccTGGCACAAAGCAACCACCAACGGATGGTCGGTGACTCGTACAGCAAGCGAATCAAAGCGGCCCGTCAACCAGCGCGGTGTTGTCGCAGGCGCGGGAAAGACAAAGGTGACAGGACCTGGCCAGCGTGAAACCATTTTTTCCCGCTggc                                                                                                                                     <  1:708478/150‑3 (MQ=255)
              cTGGCACAAAGCAACCACCAACGGATGGTCGGTGACTCGTACAGCAAGCGAATCAAAGCGGCCCGTCAACCAGCGCGGTGTTGTCGCATGCGCGGGAAAGACAAAGGTGACAGGACCTGGCCAGCGTGAAACCATTTTTTCCCGCTggcc                                                                                                                                    <  2:129993/150‑4 (MQ=255)
              cTGGCACAAAGCAACCACCAACGGATGGTCGGTGACTCGTACAGCAAGCGAATCAAAGCGGCCCGTCAACCAGCGCGGTGTTGTCGCAGGCGCGGGAAAGACAAAGGTGCCAGGACCTGGCCAGCGGGAAAAAATGGTTTCACGCTGCAc                                                                                                                                    >  2:1007180/1‑150 (MQ=255)
               tGGCACAAAGCAACCACCAACGGATGGTCGGTGACTCGTACAGCAAGCGAATCAAAGCGGCCCGTCAACCAGCGCGGTGTTGTCGCAGGCGCGGGAAAGACAAAGGTGACAGGACCTGGCCAGCGGGAAAAAATGGTTTCACGCTGCACg                                                                                                                                   <  1:1228884/150‑1 (MQ=255)
                  cacaAAGCAACCACCAACGGATGGTCGGTGACTCGTACAGCAAGCGAATCAAAGCGGCCCGTCAACCAGCGCGGTGTTGTCGCAGGCGCGGGAAAGACAAAGGTGACAGGACCTGGCCAGCGGGAAAAAATGGTTTCACGCTGCGCgtca                                                                                                                                >  1:1044359/1‑150 (MQ=255)
                  cacaAAGCAACCACCAACGGATGGTCGGTGACTCGTACAGCAAGCGAATCAAAGCGGCCCGTCAACCAGCGCGGTGTTGTCGCAGGCGCGGGAAAGACAAAGGTGACAGGACCTGGCCAGCGGGAAAAAATGGTTTCACGCTGCACgtca                                                                                                                                >  1:1402392/1‑150 (MQ=255)
                  cacaAAGCAACCACCAACGGATGGTCGGTGACTCGTACAGCAAGCGAATCAAAGCGGCCCGTCAACCAGCGCGGTGTTGTCGCAGGCGCGGGAAAGACAAAGGTGACAGGACCTGGCCAGCGGGAAAAAATGGTTTCACGCTGCACgtca                                                                                                                                >  2:978966/1‑150 (MQ=255)
                        gCAACCACCAACGGATGGTCGGTGACTCGTACAGCAAGCGAATCAAAGCGGCCCGTCAACCAGCGCGGTGTTGTCGCAGGCGCGGGAAAGACAAAGGTGACAGGACCTGGCCAGCGGGAAAAAATGGTTTCACGCTGCACGTCAGTCAAc                                                                                                                          <  1:910776/150‑1 (MQ=255)
                        gCAACCACCAACGGATGGTCGGTGACTCGTACAGCAAGCGAATCAAAGCGGCCCGTCAACCAGCGCGGTGTTGTCGCAGGCGCGGGAAAGACAAAGGTGACAGGACCTGGCCAGCGGGAAAAAATGGTTTCACGCTGCACGTCAGTCAAc                                                                                                                          <  2:977976/150‑1 (MQ=255)
                            ccaccaACGGATGGTCGGTGACTCGTACAGCAAGCGAATCAAAGCGGCCCGTCAACCAGCGCGGTGTTGTCGCAGGCGCGGGAAAGACAAAGGTGACAGGACCTGGCCAGCGGGAAAAAATGGTTTCACGCTGCACGTCAGTCAACATgg                                                                                                                      <  1:195891/150‑1 (MQ=255)
                               ccaACGGATGGTCGGTGACTCGTACAGCAAGCGAATCAAAGCGGCCCGTCAACCAGCGCGGTGTTGTCGCAGGCGCGGGAAAGACAAAGGTGACAGGACCTGGCCAGCGGGAAAAAATGGTTTCACGCTGCACGTCAGTCAACATGgtgt                                                                                                                   >  1:307869/1‑150 (MQ=255)
                                  aCGGATGGTCGGTGACTCGTACAGCAAGCGAATCAAAGCGGCCCGTCAACCAGCGCGGTGTTGTCGCAGGCGCGGGAAAGACAAAGGTGACAGGACCTGGCCAGCGGGAAAAAATGGTTTCACGCTGCACg                                                                                                                                   >  1:244809/1‑131 (MQ=255)
                                  aCGGATGGTCGGTGACTCGTACAGCAAGCGAATCAAAGCGGCCCGTCAACCAGCGCGGTGTTGTCGCAGGCGCGGGAAAGACAAAGGTGACAGGACCTGGCCAGCGGGAAAAAATGGTTTCACGCTGCACg                                                                                                                                   <  2:244809/131‑1 (MQ=255)
                                  aCGGATGGGCGGTGACTCGTACAGCAAGCGAATCAAAGCGGCCCGTCAACCAGCGCGGTGTTGTCTCAGGCGCTGGAAAGACAAAGGTGACAGGTCCTGGCCAGCGGGAAAATATGGTTTCACGCTGCACGTCAGTCAAGATGGTGtcat                                                                                                                <  2:316869/150‑1 (MQ=255)
                                    ggATGGTCGGTGACTCGTACAGCAAGCGAATCAAAGCGGCCCGTCAACCAGCGCGGTGTTGTCGCAGGCGCGGGAAAGACAAAGGTGACAGGACCTGGCCAGCGGGAAAAAATGGTTTCACGCTGCACGTCAGTCAACATGGTGtcatca                                                                                                              <  1:978966/150‑1 (MQ=255)
                                     gATGGTCGGTGACTCGTACAGCAAGCGAATCAAAGCGGCCCGTCAACCAGCGCGGTGTTGTCGCAGGCGCGGGAAAGACAAAGGTGACAGGACCTGGCCAGCGGGAAAAAATGGTTTCACGCTGCACGTCAGTCAACATGGTGTCATCaa                                                                                                             >  1:1172821/1‑150 (MQ=255)
                                           cGGTGACTCGTACAGCAAGCGAATCAAAGCGGCCCGTCAACCAGCGCGGTGTTGTCGCAGGCGCGGGAAAGACAAAGGTGACAGGACCTGGCCAGCGGGAAAAAATGGTTTCACGCTGCACGTCAGTCAACATGGTGTCATCAATATAgg                                                                                                       >  2:93749/1‑150 (MQ=255)
                                               gACTCGTACAGCAAGCGAATCAAAGCGGCCCGTCAACCAGCGCGGTGTTGTCGCAGGCGCGGGAAAGACAAAGGTGACAGGACCTGGCCAGCGGGAAAAAATGGTTTCACGCTGCACGTCAGTCAACATGGTGTCATCAATATAGGGttt                                                                                                   >  2:284246/1‑150 (MQ=255)
                                                   cGTACAGCAAGCGAATCAAAGCGGCCCGTCAACCAGCGCGGTGTTGTCGCAGGCGCGGGAAAGACAAAGGTGACAGGACCTGGCCAGTGGGAAAAAATGGTTTCACGCTGCACGTCAGTCAACATGGTGTCATCAATATAGGGTTTAAgc                                                                                               <  1:139449/150‑1 (MQ=255)
                                                   cGTACAGCAAGCGAATCAAAGCGGCCCGTCAACCAGCGCGGTGTTGTCGCAGGCGCGGGAAAGACAAAGGTGACAGGACCTGGCCAGCGGGAAAAAATGGTTTCACGCTGCACGTCAGTCAACATGGTGTCATCAATATAGGGTTTAAgc                                                                                               <  1:92142/150‑1 (MQ=255)
                                                   cGTACAGCAAGCGAATCAAAGCGGCCCGTCAACCAGCGCGGTGTTGTCGCAGGCGCGGGAAAGACAAAGGTGACAGGACCTGGCCAGCGGGAAAAAATGGTTTCACGCTGCACGTCAGTCAACATGGTGTCATCAATATAGGGTTTAAgc                                                                                               <  1:97748/150‑1 (MQ=255)
                                                     tACAGCAAGCGAATCAAAGCGGCCCGTCAAGCAGCGCGGTGTTGTCGCAGGCGCGGGAAAGACAAAGGTGACAGGACCTGGCCAGCGGGAAAAAATGGTTTCACGCTGCACGTCAGTCAACATGGTGTCATCAATATAGGGTTTAAgctg                                                                                             <  2:1055760/150‑1 (MQ=255)
                                                     tACAGCAAGCGAATCAAAGCGGCCCGTCAACCAGCGCGGTGTTGTCGCAGGCGCGGGAAAGACAAAGGTGACAGGACCTGGCCAGCGGGAAAAAATGGTTTCACGCTGCACGTCAGTCAACATGGTGTCATCAATATAGGGTTTAAgctg                                                                                             >  1:629459/1‑150 (MQ=255)
                                                      aCAGCAAGCGAATCAAAGCGGCCCGTCAACCAGCGCGGTGTTGTCGCAGGCGCGGGAAAGACAAAGGTGACAGGACCTGGCCAGCGGGTAAATATGGTTTCACGCTGCACGTCAGTCAACATGGGGTCCTCAATATAGGGTTTAAgctgc                                                                                            >  2:1389295/1‑150 (MQ=255)
                                                      aCAGCAAGCGAATCAAAGCGGCCCGTCAACCAGCGCGGTGTTGTCGCAGGCGCGGGAAAGACAAAGGTGACAGGACCTGGCCAGCGGGAAAAAATGGTTTCACGCTGCACGTCAGTCAACATGGTGTCATCAATATAGGGTTTAAgctgc                                                                                            >  2:250683/1‑150 (MQ=255)
                                                       cAGCAAGCGAATCAAAGCGGCCCGTCAACCAGCGCGGTGTTGTCGCAGGCGCGGGAAATACAAAGGTGACAGGACCTGGCCAGCGGGAAAAAATGGTTTCACGCTGCACGTCAGTCAACATGGTGTCATCAATATAGGGTTTAAgctgct                                                                                           <  2:307869/150‑1 (MQ=255)
                                                         gCAAGCGAATCAAAGCGGCCCGTCAACCAGCGCGGTGTTGTCGCAGGCGCGGGAAAGACAAAGGTGACAGGACCTGGCCAGCGGGAAAAAATGGTTTCACGCTGCACGTCAGTCAACATGGTGTCATCAATATAGGGTTTAAGCTGCTCg                                                                                         >  2:680434/1‑150 (MQ=255)
                                                         gCAAGCGAATCAAAGCGGCCCGTCAACCAGCGCGGTGTTGTCGCAGGCGCGGGAAAGACAAAGGTGACAGGACCTGGCCAGCGGGAAAAAATGGTTTCACGCTGCACGTCAGTCAACATCGTGTCATCAATATAGGGTTTAAGCTGCACg                                                                                         >  1:608166/1‑150 (MQ=255)
                                                            aGCGAATCAAAGCGGCCCGTCAACCAGCGCGCTGTTGTCGCAGGCGCGGGAAAGACAAAGGTGACAGGACCTGCCCAGCGGGAAAAAATGGTTTCACGCTGCACGTCAGTCAACATGGTGTCATCAATATAGGGTTTAAGCTGCTCCTaa                                                                                      >  2:1289820/1‑150 (MQ=255)
                                                             gCGAATCAAAGCGGCCCGTCAACCAGCGCGGTGTTGTCGCAGGCGCTGGAAAGACAAAGGTGACAGGACCTGGCCAGCGGGAAAAAATGGTTTCACGCTGCACGTCAGTGAACATGGTGTCATCAATATAGGGTTTAAGCTGCTCGTAAt                                                                                     <  2:284311/150‑1 (MQ=255)
                                                              cGAATCAAAGCGGCCCGTCAACCAGCGCGGTGTTGTCGCAGGCGCGGGAAAGACAAAGGTGACAGGACCTGGCCAGCGGGAAAAAATGGTTTCACGCTGCACGTCAGTCAACATGGTGTCATCAATATAGGGTTTAAGCTGCTCGTAAtt                                                                                    >  1:912456/1‑150 (MQ=255)
                                                                 aTCAAAGCGGCCCGTCAACCAGCGCGGTGTTGTCGCAGGCGCGGGAAAGACAAAGGTGACAGGACCTGGCCAGCGGGAAAAAATGGTTTCACGCTGCACGTCAGTCAACATGGTGTCATCAATATAGGGTTTAAGCTGCTCGTAATAtgc                                                                                 <  2:808702/150‑1 (MQ=255)
                                                                    aaaGCGGCCCGTCAACCAGCGCGGTGTTGTCGCAGGCGCGGGAAAGACAAAGGTGACAGGACCTGGCCAGCTGGAAAAAATGGTTTCACGCTGCACGTCAGTCAACATGGTGTCATCAATATAGGGTTTAAGCTGCTCGTAATTtgctgc                                                                              <  2:216688/150‑1 (MQ=255)
                                                                    aaaGCGGCCCGTCAACCAGCGCGGTGTTGTCGCAGGCGCGGGAAAGACAAAGGTGACAGGACCTGGCCAGCGGGAAAAAATGGTTTCACGCTGCACGTCAGTCAACATGGTGTCATCAATATAGGGTTTAAGCTGCTCGTAATTtgctgc                                                                              >  2:1497875/1‑150 (MQ=255)
                                                                     aaGCGGCCCGTCAACCAGCGCGGTGTTGTCGCAGGCGCGGGAAAGACAAAGGTGACAGGACCTGGCCAGCGGGAAAAAATGGTTTCACGCTGCACGTCAGTCAACATGGTGTCATCAATATAGGGTTTAAGCTGCTCGTAATTTGCTGCg                                                                             <  1:131421/150‑1 (MQ=255)
                                                                         ggCCCGTCAACCAGCGCGGTGTTGTCGCAGGCGCGGGAAAGACAAAGGTGACAGGACCTGGCCAGCGGGAAAAAATGGTTTCACGCTGCACGTCAGTCAACATGGTGTCATCAATATAGGGTTTAAGCTGCTCGTAATTTGCTGCGATTa                                                                         >  2:671336/1‑150 (MQ=255)
                                                                         ggCCCGTCAACCAGCGCGGTGTTGTCGCAGGCGCGGGAAAGACAAAGGTGACAGGACCTGGCCAGCGGGAAAAAATGGTTTCACGCTGCACGTCAGTCAACATGGTGTCATCAATATAGGGTTTAAGCTGCTCGTAATTTGCTGCGATTa                                                                         <  1:1282653/150‑1 (MQ=255)
                                                                            ccGTCAACCAGCGCGGTGTTGTCGCAGGCGCGGGAAAGACAAAGGTGACAGGACCTGGCCAGCGGGAAAAAATGGTTTCACGCTGCACGTCAGTCAACATGGTGTCATCAATATAGGGTTTAAGCTGCTCGTAATTTGCTGCGATTaaaa                                                                      >  1:416338/1‑150 (MQ=255)
                                                                                   ccAGCGCGCTGTTGTCCCAGGCGCGGGAAATACAAAGGTGACAGGACCTGGCCAGCGGGAAAAAATGGTTTCACGCTGCACGTCACTCAACATGGTGTCATCAATATAGGGTTTAAGCTGCTCGTAATTTGCTGCGATTAAAATCAGccc                                                               >  2:1129491/1‑150 (MQ=255)
                                                                                     aGCGCGGTGTTGTCGCAGGCGCGGGAAAGACAAAGGTGACAGGACCTGGCCAGCGGGAAAAAATGGTTTCACGCTGCACGTCAGTCAACATGGAGTCATCAATATAGGGTTTAAGCTGCTCGTAATTTGCTGCGATTAAAATCAGCCCCt                                                             >  1:7992/1‑150 (MQ=255)
                                                                                             gttgtCGCAGGCGCGGGAAAGACAAAGGTGACAGGACCTGGCCAGCGGGAAAAAATGGTTTCACGCTGCACGTCAGTCAACATGGTGTCATCAATATAGGGTTTAAGCTGCTCGTAATTTGCTGCGATTAAAATCAGCCCCTTATCAAcc                                                     >  2:1488482/1‑150 (MQ=255)
                                                                                              ttgtCGCAGGCGCGGGAAAGACAAAGGTGACAGGACCTGGCCAGCGGGAAAAAATGGTTTCACGCTGCACGTCAGTCAACATGGTGTCATCAATATAGGGTTTAAGCTGCTCGTAATTTGCTGCGATTAAAATCAGCCCCTTATCAACCg                                                    <  1:1497875/150‑1 (MQ=255)
                                                                                              ttgtCGCAGGCGCGGGAAAGACAAAGGTGACAGGACCTGGCCAGCGGGAAAAAATGGTTTCACGCTGCACGTCAGTCAACATGGTGTCATCAATATAGGGTTTAAGCTGCTCGTAATTTGCTGCGATTAAAATCAGCCCCTTATCAACCg                                                    <  1:680434/150‑1 (MQ=255)
                                                                                                gtCGCAGGCGCGGGAAAGACAAAGGTGACAGGACCTGGCCAGCGGGAAAAAATGGTTTCACGCTGCACGTCAGTCAACATGGTGTCATCAATATAGGGTTTAAGCTGCTCGTAATTTGCTGCGATTAAAATCAGCCCCTTATCAACCGGa                                                  >  2:518727/1‑150 (MQ=255)
                                                                                                  cGCAGGCGCGGGAAAGACAAAGGTGACAGGACCTGGCCAGCGGGAAAAAATGGTTTCACGCTGCACGTCAGTCAACATGGTGTCATCAATATAGGGTTTAAGCTGCTCGTAATTTGCTGCGATTAAAATCAGCCCCTTATCAACCg                                                    <  2:1136321/146‑1 (MQ=255)
                                                                                                  cGCAGGCGCGGGAAAGACAAAGGTGACAGGACCTGGCCAGCGGGAAAAAATGGTTTCACGCTGCACGTCAGTCAACATGGTGTCATCAATATAGGGTTTAAGCTGCTCGTAATTTGCTGCGATTAAAATCAGCCCCTTATCAACCg                                                    >  1:1136321/1‑146 (MQ=255)
                                                                                                   gCAGGCGCGGGAAAGACAAAGGTGACAGGACCTGGCCAGCGGGAAAAAATGGTTTCACGCTGCACGTCAGTCAACATGGTGTCATCAATATAGGGTTTAAGCTGCTCGTAATTTGCTGCGATTAAAATCAGCCCCTTATCAACCGGACGc                                               >  1:375965/1‑150 (MQ=255)
                                                                                                   gCAGGCGCGGGAAAGACAAAGGTGACAGGACCTGGCCAGCGGGAAAAAATGGTTTCACGCTGCACGTCAGTCAACATGGTGTCATCAATATAGGGTTTAAGCTGCTCGTAATTTGCTGCGATTAAAATCAGCCCCTTATCAACCGGACGc                                               <  1:749130/150‑1 (MQ=255)
                                                                                                      ggCGCGGGAAAGACAAAGGTGACAGGACCTGGCCAGCGGGAAAAAATGGTTTCACGCTGCACGTCAGTCAACATGGTGTCATCAATATAGGGTTTAAGCTGCTCGTAATTTGCTGCGATTAAAATCAGCCCCTTATCAACCGGACGCTGt                                            <  2:1269925/150‑1 (MQ=255)
                                                                                                        cgcgGGAAAGACAAAGGTGACAGGACCTGGCCAGCGGGAAAAAATGGTTTCACGCTGCACGTCAGTCAACATGGTGTCATCAATATAGGGTTTAAGCTGCTCGTAATTTGCTGCGATTAAAATCAGCCCCTTAt                                                          >  2:485419/1‑134 (MQ=255)
                                                                                                        cgcgGGAAAGACAAAGGTGACAGGACCTGGCCAGCGGGAAAAAATGGTTTCACGCTGCACGTCAGTCAACATGGTGTCATCAATATAGGGTTTAAGCTGCTCGTAATTTGCTGCGATTAAAATCAGCCCCTTAt                                                          <  1:485419/134‑1 (MQ=255)
                                                                                                         gcgGGAAAGACAAAGGGTACAGGACCTGGCCAGCTGGAAAAAATGGTTTCACGCTGCACGTCAGTCAACATGGTGTCATCAATATAGGGTTTAAGCTGCTCGTAATTTGCTGCGATTAAAATCAGCCCCTTATCAACCGGACGCTGtttt                                         <  2:973601/150‑1 (MQ=255)
                                                                                                          cgGGAAAGACAAAGGTGACAGGACCTGGCCAGCGGGAAAAAATGGTTTCACGCTGCACGTCAGTCAACATGGTGTCATCAATATAGGGTTTAAGCTGCTCGTAATTTGCTTCGATTAAAATCAGCCCCTTATCAACCGGACGCTGTTTTa                                        <  1:1213572/150‑1 (MQ=255)
                                                                                                            ggAAAGACAAAGGTGACAGGACCTGGCCAGCGGGAAAAAATGGTTTCACGCTGCACGTCAGTCAACATGGTGTCATCAATATAGGGTTTAAGCTGCTCGTAATTTGCTGCGATTAAAATCAGCCCCTTATCAACCGGACGCTGTTTTAAc                                      <  1:93749/150‑1 (MQ=255)
                                                                                                            ggAAAGACAAAGGTGACAGGACCTGGCCAGCGGGAAAAAATGGTTTCACGCTGCACGTCAGTCAACATGGTGTCATCAATATAGGGTTTAAGCTGCTCGTAATTTGCTGCGATTAAAATCAGCCCCTTATCAACCGGACGCTGTTTTAAc                                      <  2:1402392/150‑1 (MQ=255)
                                                                                                             gAAAGACAAAGGTGACAGGACCTGGCCAGCGGGAAAAAATGGTTTCACGCTGCACGTCAGTCAACATGGTGTCATCAATATAGGGTTTAAGCTGCTCGTAATTTGCTGCGATTAAAATCAGCCCCTTATCAACCGGACGCTGTTTTAACt                                     >  1:39972/1‑150 (MQ=255)
                                                                                                              aaaGACAAAGGTGACAGGACCTGGCCAGCGGGAAAAAATGGTTTCACGCTGCACGTCAGTCAACATGGTGTCATCAATATAGGGTTTAAGCTGCTCGTAATTTGCTGCGATTAATATCAGCCCCTTATCAACCGGACGCTGTTTTAACTc                                    <  1:1047662/150‑1 (MQ=255)
                                                                                                              aaaGACAAAGGTGACAGGACCTGGCCAGCGGGAAAAAATGGTTTCACGCTGCACGTCAGTCAACATGGTGTCATCAATATAGGGTTTAAGCTGCTCGTAATTTGCTGCGATTAAAATCAGCCCCTTATCAACCGGACGCTGTTTTAACTc                                    <  2:935372/150‑1 (MQ=255)
                                                                                                               aaGACAAAGGTGACAGGACCTGGCCAGCGGGAAAAAATGGTTTCACGCTGCACGTCAGTCAACATGGTGTCATCAATATAGGGTTTAAGCTGCTCGTAATTTGCTGCGATTAAAATCAGCCCCTTATCAACCGGACGCTGTTTTAACTca                                   >  2:489049/1‑149 (MQ=255)
                                                                                                                 gACAAAGGTGACAGGACCTGGCCAGCGGGAAAAAATGGTTTCACGCTGCACGTCAGTCAACATGGTGTCATCAATATAGGGTTTAAGCTGCTCGTAATTTGCTGCGATTAAAATCAGCCCCTTATCAACCGGACGCTGTTTTAACTCCaa                                 >  2:1027445/1‑150 (MQ=255)
                                                                                                                  aCAAAGGTGACAGGACCTGGCCAGCGGGAAAAAATGGTTTCACGCTGCACGTCAGTCAACATGGTGTCATCAATATAGGGTTTAAGCTGCTCGTAATTTGCTGCGATTAAAATCAGCCCCTTATCAACCGGACGCTGTTTTAACTCCAAc                                <  2:994437/150‑1 (MQ=255)
                                                                                                                   cAAAGGTGACAGGACCTGGCCAGCGGGAAAAAATGGTTTCACGCTGCACGTCAGTCAACATGGTGTCATCAATATAGGGTTTAAGCTGCTCGTAATTTGCTGCGATTAAAATCAGCCCCTTATCAACCGGACGCTGTTTTAACTCCAACa                               >  2:733478/1‑150 (MQ=255)
                                                                                                                   cAAAGGTGACAGGACCTGGCCAGCGGGAAAAAATGGTTTCACGCTGCACGTCAGTCAACATGGTGTCATCAATATAGGGTTTAAGCTGCTCGTAATTTGCTGCGATTAAAATCAGCCCCTTATCAACCGGACGCTGTTTTAACTCCAACa                               >  1:1116925/1‑150 (MQ=255)
                                                                                                                         tGACAGGACCTGGCCAGCGGGAAAAAATGGTTTCACGCTGCACGTCAGTCAACATGGTGTCATCAATATAGGGTTTAAGCTGCTCGTAATTTGCTGCGATTAAAATCAGCCCCTTATCAACCGGACGCTGTTTTAACTCCAACAGTCGCa                         <  2:1005871/150‑1 (MQ=255)
                                                                                                                         tGACAGGACCTGGCCAGCGGGAAAAAATGGTTTCACGCTGCACGTCAGTCAACATGGTGTCATCAATATAGGGTTTAAGCTGCTCGTAATTTGCTGCGATTAAAATCAGCCCCTTATCAACCGGACGCTGTTTTAACTCCAACAGTCGCa                         <  1:875738/150‑1 (MQ=255)
                                                                                                                          gACAGGACCTGGCCAGCGGGAAAAAATGGTTTCACGCTGCACGTCAGTCAACATGGTGTCATCAATATAGGGTTTAAGCTGCTCGTAATTTGCTGCGATTAAAATCAGCCCCTTATCAACCGGACGCTGTTTTAACTCCAACAGTCGCAt                        <  1:613589/150‑1 (MQ=255)
                                                                                                                          gACAGGACCTGGCCAGCGGGAAAAAATGGTTTCACGCTGCACGTCAGTCAACATGGTGTCATCAATATAGGGTTTAAGCTGCTCGTAATTTGCTGCGATTAAAATCAGCCCCTTATCAACCGGACGCTGTTTTAACTCCAACAGTCGCAt                        >  2:805146/1‑150 (MQ=255)
                                                                                                                             aGGACCTGGCCAGCGGGAAAAAATGGTTTCACGCTGCACGTCAGTCAACATGGTGTCATCAATATAGGGTTTAAGCTGCTCGTAATTTGCTGCGATTAAAATCAGCCCCTTATCAACCGGACGCTGTTTTAACTCCAACAGTCGCATCAc                     >  1:1134486/1‑150 (MQ=255)
                                                                                                                              ggACCTGGCCAGCGGGAAAAAATGGTTTCACGCTGCACGTCAGTCAACATGGTGTCATCAATATAGGGTTTAAGCTGCTCGTAATTTGCTGCGATTAAAATCAGCGCCTTATCAACCGGACGCTGTTTTCACTCCGACAGTCGCATCagt                    >  2:1053930/1‑148 (MQ=255)
                                                                                                                              ggACCTGGCCAGCGGGAAAAAATGGTTTCACGCTGCACGTCAGTCAACATGGTGTCATCAATATAGGGTTTAAGCTGCTCGTAATTTGCTGCGATTAAAATCAGCCCCTTATCAACCGGACGCTGTTTTAACTCCAACAGTCGCATCAct                    >  2:313982/1‑150 (MQ=255)
                                                                                                                                 ccTGGCCAGCGGGAAAAAATGGTTTCACGCTGCACGTCAGTCAACATGGTGTCATCAATATAGGGTTTAAGCTGCTCGTAATTTGCTGCGATTAAAATCAGCCCCTTATCAACCGGACGCTGTTTTAACTCCAACAGTCGCATCActgct                 >  2:1526230/1‑150 (MQ=255)
                                                                                                                                  cTGGCCAGCGGGAAAAAATGGTTTCACGCTGCACGTCAGTCAACATGGTGTCATCAATATAGGGTTTAAGCTGCTCGTAATTTGCTGCGATTAAAATCAGCCCCTTATCAACCGGACGCTGTTTTAACTCCAACAGTCGCATCActgctg                >  1:584741/1‑150 (MQ=255)
                                                                                                                                   tGGCCAGCGGGAAAAAATGGTTTCACGCTGCACGTCAGTCAACATGGTGTCATCAATATTGGGTTTAAGCTGCTCGTAATTTGCTGCGATTAAAATCAGCCCCTTATCAACCGGACGCTGTTTTAACTCCAACAGTCGCATCACTGCTGt               >  2:1350968/1‑150 (MQ=255)
                                                                                                                                   cgtgcAGCGTGAAACCATTTTTTCCCGCTGCACGTCAGTCAACATGGTGTCATCAATATAGGGTTTAAGCTGCTCGTAATTTGCTGCGATTAAAATCAGCCCCTTATCAACCGGACGCTGTTTTAACTCCAACAGTCGCATCACTGCTGt               >  1:775977/5‑150 (MQ=255)
                                                                                                                                    ggCCAGCGGGAAAAAATGGTTTCACGCTGCACGTCAGTCAACATGGTGTCATCAATATAGGGTTTAAGCTGCTCGTAATTTGCTGCGATTAAAATCAGCCCCTTATCAACCGGACGCTGTTTTAACTCCAACAGTCGCATCACTGCTGtt              >  1:516254/1‑150 (MQ=255)
                                                                                                                                     tgcAGCGTGAAACCATTTTTTCCCGCTGCACGTCAGTCAACATGGTGTCATCAATATAGGGTTTAAGCTGCTCGTAATTTGCTGCGATTAAAATCAGCCCCTTATCAACCGGACGCTGTTTTAACTCCAACAGTCGCATCACTGCTGttt             >  1:549029/3‑150 (MQ=255)
                                                                                                                                      ccAGCGGGAAAAAATGGTTTCACGCTGCACGTCAGTCAACATGGTGTCATCAATATAGGGTTTAAGCTGCTCGTAATTTGCTGCGATTAAAATCAGCCCCTTATCAACCGGACGCTGTTTTAACTCCAACAGTCGCATCACTGCTGTTTc            >  2:944783/1‑150 (MQ=255)
                                                                                                                                       cAGCGGGAAAAAATGGTTTCACGCTGCACGTCAGTCAACATGGTGTCATCAATATAGGGTTTAAGCTGCTCGTAATTTGCTGCGATTAAAATCAGCCCCTTATCAACCGGACGCTGTTTTAACTCCAACAGTCGCATCACTGCTGTTTCg           <  2:908710/150‑1 (MQ=255)
                                                                                                                                       cAGCGGGAAAAAATGGTTTCACGCTGCACGTCAGTCAACATGGTGTCATCAATATAGGGTTTAAGCTGCTCGTAATTTGCTGCGATTAAAATCAGCCCCTTATCAACCGGACGCTGTTTTAACTCCAACAGTCGCATCACTGCTGTTTCg           <  2:327755/150‑1 (MQ=255)
                                                                                                                                            ggAAAAAATGGTTTCACGCTGCACGTCAGTCAACATGGTGTCATCAATATAGGGTTTAAGCTGCTCGTAATTTGCTGCGATTAAAATCAGCCCCTTATCAACCGGACGCTGTTTTAAc                                      >  1:1312290/1‑118 (MQ=255)
                                                                                                                                            ggAAAAAATGGTTTCACGCTGCACGTCAGTCAACATGGTGTCATCAATATAGGGTTTAAGCTGCTCGTAATTTGCTGCGATTAAAATCAGCCCCTTATCAACCGGACGCTGTTTTAAc                                      <  2:1312290/118‑1 (MQ=255)
                                                                                                                                            ggAAAAAATGGTTTCACGCTGCACGTCAGTCAACATGGTGTCATCAATATAGGGTTTAAGCTGCTCGTAATTTGCTGCGATTAAAATCAGCCCCTTATCAACCGGACGCTGTTTTAACTCCAACAGTCGCATCACTGCTGTTTCGCTATc      >  2:311791/1‑150 (MQ=255)
                                                                                                                                               aaaaaTGGTTTCACGCTGCACGTCAGTCAACATGGTGTCATCAATATAGGGTTTAAGCTGCTCGTAATTTGCTGCGATTAAAATCAGCCCCTTATCAACCGGACGCGGTTTTAACTCCAACAGTCGCATCACTGCTGTTTCGCTATCAgg   >  2:553810/1‑150 (MQ=255)
                                                                                                                                                aaaaTGGTTTCACGCTGCACGTCAGTCAACATGGTGTCATCAATATAGGGTTTAAGCTGCTCGTAATTTGCTGCGATTAAAATCAGCCCCTTATCAACCGGACGCTGTTTTAACTCCAACAGTCGCATCACTGCTGTTTCGCTATCAGGa  >  2:354057/1‑150 (MQ=255)
                                                                                                                                                 |                                                                                                                                                    
GGTTTACCATAAGCCTGGCACAAAGCAACCACCAACGGATGGTCGGTGACTCGTACAGCAAGCGAATCAAAGCGGCCCGTCAACCAGCGCGGTGTTGTCGCAGGCGCGGGAAAGACAAAGGTGACAGGACCTGGCCAGCGGGAAAAAATGGTTTCACGCTGCACGTCAGTCAACATGGTGTCATCAATATAGGGTTTAAGCTGCTCGTAATTTGCTGCGATTAAAATCAGCCCCTTATCAACCGGACGCTGTTTTAACTCCAACAGTCGCATCACTGCTGTTTCGCTATCAGGA  >  NC_000913/3431015‑3431308

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